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Crack  length,  (in  or  cm) 

Initial  crack  length,  (in  or  cm) 

Horizontal  line;  eutectoid  reaction,  (°F  or  °C) 

Critical  crack  length,  (in  or  cm) 

Final  crack  length  (precracking),  (in  or  cm) 

Crack  length  ratio,  or  dimensionless  crack  length, 
(in/in,  or  cm/cm) 

Thickness  of  specimen,  (in  or  cm) 

Slope  of  regression  line 

Compliance  or  reciprocal  of  stiffness,  (in/lb,  or  cm/nt) 
Constant  in  Paris  equation  (empirical) 

Total  crack  length,  (in  or  cm) 

Dimensionless  compliance 

Crack  opening  displacement,  also  6,  (in  or  cm) 

Compact  tension 

Crack  Growth  Rate,  also  H,  (in/cycle,  or  mm/cycle) 

2 2 

Young's  modulus,  (Ib/in  , or  nt/cm  ) 

Weibull  threshold  parameter 
Elastic  elongation,  (in  or  cm) 

Dimensionless  compliance,  also  CLB 
/.K/Kb, (ksi/ksi , or  MPa/MPa) 

K calibration  factor,  (in/in,  or  cm/cm) 
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f Function 

G,  Crack  extension  force,  (lb/in,  or  nt/cm),  also  energy 

release  rate  (Ib-in/in?,  or  nt-cm/cnr)  normal  to  flaw  pro- 
pagation 

Gc  Critical  crack  extension  force,  (ksi  or  MPa) 

H Oa/ON,  (in/cycle,  or  mn/cycle) 

H/W  CT  specimen  heighth  to  width  ratio,  (in/in,  or  cm/cm) 

K Stress  intensity,  (ksi  /Tn,  or  MPa  /m) 

K.  Final  stress  intensity  for  fitting  relations,  (ksi  /Tn, 

or  MPa  /in) 

K Critical  plane  stress  fracture  toughness,  (ksi  /Tn,  or 

MPa  /m) 

Kj  Plane  strain  stress  intensity,  (ksi  /Tn,  or  MPa  /in) 

K.  Critical  plane  strain  fracture  toughness,  (ksi  /ill,  or 

lc  MPa  /T) 

K u Maximum  stress  intensity,  (ksi  /Tn,  or  MPa  *11) 

max 

K . Minimum  stress  intensity,  (ksi  /Tn,  or  MPa  /m) 
min 

«th  Stress  intensity  threshold,  (ksi  /Tn,  or  MPa  »/in) 

k Wei  bull  parameter  for  curve  shape 

ksi  Kilo-pounds  per  square  inch  (1  ksi  = 6.895  MPa) 

LVDT  Linear  variable  differential  transformer 

MPa  Mega  Pascal,  or  106  Pascal  (unit  of  pressure  or  stress) 

M Spring  constant,  (in/lb,  or  cm/nt) 

m Fmpirical  Rate  exponents 

P Load,  (lb  or  nt) 


f 


R Stress  ratio  (Sfnin/Sniax) , (lb/in2,  or  nt/crn2) 

r Sample  coefficient  of  correlation,  also  R 

r*  Risk  function 

2 2 

5 x Stress  maximum,  (lb/in  , or  nt/cm  ) 

2 2 

^min  Stress  minimum,  (lb/in  , or  nt/cm  ) 

U Stored  linear  strain  energy,  (in-lb,  or  cm-nt) 

v Weibull  characteristic  value  for  specimen  type 

w Weibull  parameter,  1/k 

WUL  Wedge  open  loading 

Y Polynomial  expression 

A Difference; 

Range ; 2 2 

Area  of  triangle,  (in*  or  cm) 

6 COD,  deflection  of  crack  opening,  (in  or  cm) 
Elliptical  parameter,  non-dimensional 

t Strain  and  increments,  (in/in,  or  cm/cm) 

i.  True  fracture  strain,  (in/in,  or  cm/cm) 

t.  Critical  strain,  (in/in,  or  cm/cm) 

tp  Plastic  strain,  (in/in,  or  cm/cm) 

i Surface  energy,  (lb/in,  or  dyne/cm)'. 

p Radius  of  curvature,  (in  or  cm) 

w Micrometer  ( = 10"^  nt) 

v Poisson's  ratio 

Statistical  scatter  in  the  vertical  direction 
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2 2 

Uniformly  applied  normal  stress,  (lb/in  , or  nt/cm  ) 
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2 2 

True  fracture  stress,  (Ib/in  , or  nt/cm  ) 

2 2 

Critical  value  of  stress,  (lb/in  , or  nt/cm  ) 

2 2 

Tensile  strength, (lb/ in  , or  nt/cm  ) 

2 2 
Yield  strength,  (lb/in  , or  nt/cm  ) 

2 2 
Shear  stress,  (Ib/in  , or  nt/cm  ) 
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ABSTRACT 


Consumers  are  recognizing  the  need  for  close  control  in  high 
performance  designs,  materials,  and  structures  when  operating 
within  the  realm  of  rational  continuum  or  sharp  flaw  mechanics. 

This  is  evidenced  by  the  increasing  use  of  Structural  Integrity 
Technology  in  advanced  high  performance  designs.  In  order  to 
implement  a Structural  Integrity  plan,  fundamental  understanding  of 
subcritical  crack  growth  rates  is  needed.  To  achieve  this,  it  is 
essential  to  increase  the  data  base  accuracy  and  relate  physical 
boundary  condi tiojis  of  the  fatigue-crack  growth  process  to  life 
prediction  models.  The  scientific  hypothesis  aAd-li^Lj^toncJusion  - 
of  this  work  is  that  microstructure  when  acted  upon  by  load  and 
environment  determines  the  fatigue-crack  growth  rates  or  life  of  a 
component. 

The  experimental  approach  varied  the  heat  treatment  of  two 
precipitation  hardening  martensitic,  alloys,  17-4  PH  and  15-S  PH. 
Fatigue-crack  growth  data  was  correlated  between  materials  cycled 
under  wave  forms  of  10  Hz  sine  frequency  and  trapezoidal  one  minute 
hold-time  for  different  stress  range  ratios  and  environments.^' 

The  data  are  represented  through  four  parameters  of  a unique 
functional  curve,  not  a straight  line,  that  precisely  applies 
mathematical  boundary  conditions  to  the  physics  of  the  fatigue- 
crack  growth  process.  This  contemporary  Weibul 1-type  mathematical 
uethod  has  been  demonstrated  to  be  applicable  to  the  quantification 
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of  microstructural  effects  on  fatigue-crack  growth. 

It  is  the  judgnent  of  this  author  that  a key  to  future  advance- 
ment of  technology  utilization  of  fatigue-crack  growth  data  lies  in 
its  application  and  interpretation  for  structures  in  accordance  with 
physical  boundary  conditions  under  the  operational  conditions  of 
interest. 
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INTRODUCTION 

In  recent  years,  materials'  properties  and  structural  design 
requirements  have  increased  as  high  performance  technology  was  ex- 
plored and  pursued.  Using  performance  as  a standard,  materials  are 
being  used  at  their  limits  not  only  in  strength  and  toughness,  but 
also  in  fatigue  resistance.  Attempts  are  being  made  to  extend  these 
limits  through  improved  alloying,  microstructural  control,  and 
quality  control.  In  certain  areas,  changing  materials  have  given 
v/ay  to  a more  refined  approach  that  studies  component  or  structural 
operation,  expected  life,  and  maintenance  schedule  versus  its  desired 
material  response. 

These  trends  have  placed  extensive  demands  upon  the  designer  and 
materials  engineer,  for  he  must  know  with  precision  and  confidence 
the  factors  that  describe  materials,  the  operational  environments, 
and  the  anticipated  structural  response.  Achieving  this  demand  does 
not  lie  solely  in  more  thorough  data  collection  techniques,  in  data 
reduction  processes,  or  in  approaches  that  advocate  only  refinement 
of  previous  techniques  but  in  recognizing  the  physics  of  the  event. 

The  reason  is  that  neither  materials'  properties  nor  operational  en- 
vironments are  presently  subject  to  unique  description.  These  factors 
must  be  quantized  as  a function  of  the  physical  interactions  limiting 
structural  life.  The  import  of  this  approach  can  be  considered  pro- 
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portional  to  the  structure's  mission,  risk,  and  cost. 

As  a consequence  of  the  probabilistic  nature  of  materials  and 
operational  environments,  a new  design  philosophy  has  been  cultivated 
Through  a natural  progression  of  events,  (e.g.  structural  in- 
operability, failures,  and  performance  derating)  a limited  number 
of  people  have  become  advocates  of  a philosophy  termed  Structural 
Integrity  to  produce  a needed  engineering  capability  for  extending 
or  assuring  design  life.  The  uniqueness  of  this  approach  lies  in 
combining  established  disciplines  such  as  statistics,  stress  analysis 
material  science  and  fracture  mechanics,  to  produce  this  capability. 
The  combination  will  benefit  certain  disciplines  more  than  others  de- 
pending upon  the  situation,  however,  the  resultant  capability  of 
Structural  Integrity  is  the  main  concern. 

There  are  several  objectives  of  this  work  which  evaluates  the 
interaction  of  metallurgy,  stress  range,  flaw  type,  load  profile, 
and  environment  under  fatigue-crack  propagation  testing  of  precipi- 
tation hardened  (PM)  stainless  steels.  Testing  these  alloys  under 
various  conditions,  may  provide  immediate  application  of  Structural 
Integrity  to  the  U.S.  Navy  high  performance  hydrofoil  program,  PHM, 
where  major  structural  elements  are  fabricated  from  these  materials 
(Bauman,  1974). 


Chapter  2 


OBJECTIVE 

This  study  was  directed  toward  achieving  the  following  specific 
technical  objectives: 

i.  To  evaluate  the  effect  of  fatigue-crack  growth  behavior 
of  two  precipitation  hardened  (PH)  stainless  steels  (i.e. 

17-4  and  15-5)  in  two  different  heat  treated  conditions 
(i.e.  H 1050  and  H 1100  conditions.) 

ii.  To  explore  the  effect  of  wave  form  in  terms  of  increasing 
or  decreasing  the  rate  of  fatigue-crack  growth. 

Prime  importance  will  be  placed  on  the  hold-time  at  maxi- 
mum load. 

iii.  To  determine  environmental  effects  when  combined  with  wave 
form  under  hold-time  conditions. 

iv.  To  consider  the  effects  of  stress  range  ratio  with  specific 
regard  to  internal  residual  stresses, 

v.  To  provide  an  approach  and  data  to  observe  if  traditional 
methods  for  fatigue-crack  growth  testing  that  have  con- 
sidered the  influence  of  operational  load  time  profiles  are 
capable  of  describing  component  life,  maintenance,  or  mor- 
tality time  frame. 

As  indicated,  the  primary  materials  in  this  study  are  17-4  PH 
and  15-5  PH  stainless  steels.  They  are  used  in  a high  strength 
salt  water  application.  Both  base  materials  will  be  studied  over  a 
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narrow  yield  strength  level  as  controlled  by  the  final  heat  treat- 
ments. The  program  was  formulated  around  the  concepts  of  fracture 
mechanics  using  a modified  compact  tension  specimen  called  a wedge 
open  loading  (WOL)  test  specimen  of  type  1 design  (i.e.  load  normal 
to  crack  growth.) 

It  is  anticipated  that  through  this  program  the  following 
achievements  can  be  realized: 

(a) .  Provide  a satisfactory  technical  rationale  for  critical 

resolution  of  the  effect  of  environment,  wave  form,  stress 
ratio,  and  material  condition  on  fatigue-crack  growth  of 
17-4  PH  and  15-5  PH  stainless  steels. 

(b) .  Close  control  and  observation  of  microstructures,  multi- 

laboratory data  acquisition,  and  a contemporary  procedure 
for  fatigue-crack  growth  data  reduction  may  provide  a more 
decisive  mechanism  in  describing  and  relating  microstruc- 
tural  response  to  service  applications. 

(c) .  Determine  the  feasi bi 1 i ty .through  physically  definable  data 

reduction  and  fractography , of  a technique  for  assessing  a 
metals'  ability  for  stress  relaxation  in  fatigue-crack  growth 
tests  or  field  operational  applications. 
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TECHNICAL  DISCUSSION 

3.1  HISTORICAL  REVIEW  OF  FRACTURE  MECHANICS 

Considering  fracture  mechanics  as  a function  of  crack  behavior, 
metal  damage  can  be  partitioned  into  three  selected  areas  (Hoeppner, 
Fall  1975): 

1.  initiation 

2.  acceleration 

3.  unstable  crack  growth 

Traditionally,  each  of  these  areas  has  been  studied  independently 
even  though  all  areas  may  be  present  during  a single  physical 
event.  Depending  upon  the  situation,  fracture  initiation  may  carry 
the  main  import  of  crack  growth.  Assuming  this  is  appropriate,  the 
first  outgrowth  could  represent  total  failure,  however,  in  other 
situations  this  is  not  the  case.  A succinct  overview  of  these  areas 
is  warranted  considering  their  interpretation  and  interaction. 

3.1.1  Flaw  Initiation 

Usually  fracturing  begins  from  a small  discontinuity  in  what  is 
considered  a homogeneous,  continuous  material.  The  imperfection, 
termed  defect,  flaw,  inclusion,  or  void,  can  be  found  throughout  a 
structure  and  will  initiate  fracture  if  a complex  combination  of 
factors,  not  as  yet  fully  understood,  exist.  These  factors  include 
size  and  shape  of  the  imperfection,  local  stress  and  strain  fields, 
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micro-  and  macroscopic  material  characteristics,  component  load 
history,  and  environment.  These  factors  may  be  evaluated  for  the 
simple  case  but  one  must  consider  the  controlling  element,  fracture 
criterion.  This  criterion  is  the  road  map  for  overall  consideration 
of  individual  factors  that  predicate  when  fracture  will  occur. 

For  a simple  case,  the  most  widely  known  criterion  is  that  put  forth 
by  A.  A.  Griffith  (1921).  His  realization  that  strain  energy  in  a 
cracked  brittle  plate  decreases  as  the  crack  grows  has  promulgated 
mst  of  the  modem  work  in  fracture.  T he  resulting  relationship  be- 
tween increased  energy  associated  with  the  generation  of  new  crack 
surfaces  and  decrease  of  input  strain  energy  has  been  developed 
into  functional  form  and  applied  to  complex  situations  that  can  pre- 
dict final  metal  failure  beyond  that  of  basiu  brittle  materials. 

Aware  that  a.i  energy  transfer  and  balance  exists,  the  factors  af- 
fecting tne  rate  of  crack  development  or  crack  acceleration  must  be 
considered. 

2.1.2  Acceleration 

As  a crack  begins  to  grow  through  the  material  its  growth  rate 
becomes  important.  If  the  controlling  microstructural  parameters 
provide  resistance  to  crack  growth,  the  crack  acceleration  may  be  low 
enough  so  that  operational  equipment  could  be  maintained  over  its 
design  lifetime.  However,  if  crack  extension  occurs  by  bursting 
through  preferred  microstructure  paths  in  the  material,  then  the 
potential  decrease  in  component  life  imrt  he  cetr.i  dered.  Ilow  these 
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accelerations  interface  in  fatigue  loading  and  crack  growth  rates 
requires  knowledge  of  events  after  initiation  of  fracture  (Hoeppner, 
1973-1976)  and  recognition  of  impact  to  life  prediction  and  time  to 
terminal  or  unstable  fracture. 

3.1.3  Unstable  Crack  Growth 

As  crack  growth  instability  occurs  there  is  an  acceleration  to 
a constant  velocity  that  is  controlled  by  the  material's  inherent 
ability  to  absorb  strain.  In  materials  with  a high  propensity  toward 
cleavage,  a less  energetic  fracture  process,  the  velocity  is  approxi- 
mately equal  to  one-third  of  the  dilatational  wave  velocity 
(5000  ft/sec  or  1524  m/sec)  (Swedlow,  1922).  The  importance  of  the 
speed  of  unstable  crack  propagation  can  be  seen  in  the  difference 
between  static  and  dynamic  (steady  state)  stress  fields  as  applied 
to  "crack  arrestor"  design  in  wings,  fuselages,  tanks,  and  piping 
systems. 

3.2  FRACTURE  CRITERION 

As  a composite  of  loads  are  placed  on  a specimen,  a component 
with  a small  flaw,  for  example,  fracture  may  begin.  Assuming  it  has 
initiated  and  a definition  of  initiation  has  been  agreed  upon, 
Griffith  considered  two  questions; 

1)  What  type  of  mathematical  model  of  the  event  can 
be  appl ied? 

2)  What  defines  the  critical  load  that  made  the 
fracture  begin  for  a simple  model  material? 
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3.2. 1 Gri f f i th  Model 

Griffith's  mathematical  model  was  required  to  describe  a physi- 
cal situation  by  way  of  providing  necessary  information  and  be  in  a 
form  such  that  analytical  techniques  could  be  applied  to  the  model. 
Besides  geometry  and  loading  parameters,  other  factors  that  deserved 
consideration  were  material  characteristics  such  as  metallurgy  versus 
e’astic  values,  residual  stress,  load  history,  and  environment 
(Hoeppner  et  al,  i974).  However,  the  abundance  of  physical  factors 
met  with  limited  analytical  techniques  such  as  dislocation  theory, 
plasticity,  elasticity,  finite  elements  and  statistics.  Even  though 
elasticity  does  not  consider  chemical  composition,  heat  treat  or 
environment,  the  continuum  viewpoint  which  was  most  familiar  was  to 
be  the  basis  of  the  model  in  spite  of  idealization. 

The  Griffith  model  was  a flat  plate  of  linearly  elastic,  iso- 
tropic, homogeneous  material  that  followed  the  relationships  of 
classical  elastomechanics.  An  infinitely  long  and  wide  plate  con- 
taining a crack  Za  long  through  Us  thickness  was  selected.  Its 
stress  and  strain  fields  were  considered  invariant  with  respect  to 
thickness  constraint.  The  perimeter  of  the  plate  was  equipol lently 
stressed. 

3.2.2  Interpretation 

Since  the  flaw  has  been  represented  by  a perfectly  sharp  crack, 
fracturing  was  considered  a linear  extension  of  the  crack,  i.e.  2a 
to  2(a  + ''a).  The  condition  for  extension  could  have  been  a value  of 
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stress  as  identified  by  maximum  shear  or  octahedral  shear  or  some 
limit  of  strain.  However,  these  parameters  were  not  applicable  due 
to  the  stress  field's  position  with  respect  to  the  crack  tip 
(Swedlow,  1972). 

In  lieu  of  setting  governing  parameters  in  a spurious  manner  a 
key  physical  factor  and  necessary  condition  was  conservation  of 
energy.  The  function,  as  previously  stated,  maintains  a balance  of 
relative  energy  change  between  new  surface  energy  and  strain  energy 
for  a crack  incremental  extension.  As  applied  to  this  model  the  con- 
dition has  become  known  as  Griffith's  criterion  for  brittle  frac- 
ture (Griffith,  1921). 

Griffith's  work  is  based  on  three  salient  points,  definition 
of  stress  field  at  the  crack,  mode  of  crack  growth,  and  criterion. 

The  first  element  is  extracted  from  the  work  of  Inglis  (1913)  where 
he  computed  the  stresses  in  a plate  with  an  elliptically  shaped 
hole.  The  describing  equation  with  foci  of  the  ellipse  on  the 
x-axis  at  x •+  a is: 

*2.  + —XL ,-=l  (1) 

(a  cosh  i ) (a  sinh  ;'.o) 

where  t,  is  a non-negative  parameter  that  when  approaching  zero, 
shrinks  the  ellipse  to  a line  length  2a  along  the  x-axis.  Succinctly, 
Inglis  found  the  stress  field  for  various  loadings;  he  located  the 
highest  stress  at  the  end  of  the  major  x-axis.  The  stress  was  approxi- 
mated by  applied  stress  o times  2/a7i7,  where  ,■  the  radius  of  curva- 
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ture  is  expressed  as 


p = a sinh  tanh  CQ 


(2) 


For  a crack,  as  £0  goes  to  zero  all  components  of  stress  have 
nagnitudes  inversely  proportional  to  the  square  root  of  the  distance 
from  the  crack  tip.  This  characteristic  stress  is  either  infinite  or 
singular  and  the  strain  is  similar  as  calculated  by  linear  elastic- 
ity. At  the  crack  tip,  the  displacements  are  bounded  (Williams, 

1957).  Thus,  limiting  factors  of  stress  or  strain  as  fracture  cri- 
teria were  inapp'Opri ate  for  this  model. 

The  potential  strain  energy  due  to  a crack  is  mathematical ly  in- 

2 2 

tegrable,  thus  a finite  quantity  directly  proportional  to  (a  o /l), 
per  unit  plate  thickness.  Computing  the  total  potential  energy  of 
the  cracked  plate  Griffith  added  a factor  4ay  where  \ equals  surface 
energy.  Assuming  fracture  was  attributed  to  crack  extension  the  total 
energy  was  to  remain  unchanged.  The  result  is  a critical  value  of  a 
designated  by  o*  and  given  by: 


u*  = 


na  (1- 


TT 

v ) 


1/2 


for  plane  strain 


(3) 


ztL 

ii  a 


1/2 


(4) 


for  plane  stress 


Griffith's  criterion  predicts  a relationship  between  critical 
stress,  geometry,  and  elastic  properties  of  an  idealized  plate.  It 
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can  be  stated  that  a*/a~  is  equal  to  a ma  , constant. 

As  Griffith's  criterion  was  developed  and  derived  from  work 
with  glass  as  a model  material  it  was  promoted  later  for  application 
to  more  common  engineering  materials.  A discrepancy  was  observed 
when  under  restrictive  conditions  the  value  o*/a  was  not  close  to  a 
constant.  It  was  hypothesized  that  creation  of  new  surfaces  was  not 
a total  representation  of  energy  interchange  for  elastic  strain 
energy  an  aspect  Griffith  recognized  as  a physical  reality  that 
would  have  to  be  considered.  Traditionally,  Orowan  (1952)  and 
Irwin  (1948)  are  credited  with  amending  Griffith's  original  formula 
to  reflect  plasticity  conditions.  Their  modification  attempts  to 
account  for  plastic  strain  energy  and  flow  at  the  crack  tip  absorb- 
ing or  dissipating  some  potential  energy.  The  amended  formula  can 
be  stated  as: 

^aO-v2) 

where  Gjc  (1  is  the  opening  mode  normal  to  crack)  is  an  energy  re- 
lease rate  with  crack  extension,  replacing  2)  and  measured  directly 
from  a fracture  test.  Presently,  there  are  approximate  solutions  for 
the  energy  dissipation  due  to  crack  tip  plasticity,  but  because  all 
microstructural  variables  and  physical  fracture  processes  are  not 
uniquely  described,  exact  solutions  considering  material  compati- 
bility and  equilibrium  conditions  are  not  available  (Hoeppner, 
1973-1976). 
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3.2.3  Stress  Intensity  Consideration 

Basic  equations  of  linear  elastic  fracture  nechanics  relate 
Gjc  to  the  elastic  stress  and  flaw  size  conditions  necessary 
for  fracture.  These  equations  are  utilized  in  the  practical  sense 
only  for  the  case  of  plane  strain,  because  of  strong  dependence  of 
Gjc  on  specimen  dimensions  and  materials'  properties.  For  plane 
strain,  as  noted  by  Irwin  and  others,  G is  related  to  the  criti- 
cal  stress  intensity  factor  Kjc  by  K.c  = GjcE  where  E is  Young's 
Modulus  (Judy,  1972).  Using  a simple  model  from  previous  publica- 
tions of  Inglis  (1913)  and  Williams  (1957),  and  deleting  crack-tip 
plasticity  effects,  Irwin  was  able  to  analyze  the  stress  state  at 
a point  remote  from  the  crack  tip  (Figure  1).  In  these  relations, 

Kj  is  proportional  to  o/a;  with  all  other  components  of  tiie  equations 
pertaining  to  geometry.  It  is  important  to  notice  that  the  stress 
intensity  factors  are  not  dependent  on  the  coordinates,  r and  0; 
hence  they  control  the  intensity  of  the  stress  fields  and  reflect 
the  redistribution  of  stress  in  a body  due  to  a crack.  Thus,  Kj  is 
a measure  of  Intensity  of  stress  singularity,  deduced  from  an  appro- 
priate analysis  of  the  local  elastic  stress  fields  (Williams,  1957), 
and  it  has  nothing  whatever  to  do  with  fracture  of  materials  per  se. 
Simply  stated,  the  Gri ffith- Irwin  work  defines  an  equation  relating 
stress  intensity,  K,  flaw  size,  a,  and  applied  stress,  !. : 

K ; S * a a C (o ) 


Figure  1 Model  and  Equations  for  Describing  the  Stress  Field  near  the 
Crack  Tip  in  a Body.  In  these  equations  the  relation  K o/a 
is  obviQus;  all  other  parts  of  the  equation  pertain  to 
geometry. 
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where  C considers  geometry,  plasticity  approximations  and  other 
effects  (Hoeppner,  Fall  1975). 

However,  critical  values  of  K describe  characteristics  of  ma- 
terial. For  brittle  materials,  the  point  in  the  energy  balance  sys- 
tem, proposed  by  Griffith,  where  the  onset  of  crack  extension  is 
reached  is  a point,  Gj  . The  same  point,  described  by  critical  stress 
intensity,  Kj^,  has  specific  (leaning  in  that  it  is  a material  property 
measuring  propensity  of  a specific  material  to  fracture.  It  is  con- 
sidered an  invariant  quantity  within  reasonable  limits.  Similarly, 

Gjc  is  a material  property.  However,  applied  G and  K connote  totally 
different  ideas.  The  term  that  describes  available  energy  is  G where 
the  term  describing  stress  state  is  K.  The  two  terms  are  then  related 
only  at  points  where  crack  extension  takes  place  and  both  apply  only 
for  initiation  of  cracking. 

3.2.4  Fracture  Mechanics  to  Fatigue 

It  is  beneficial  to  summarize  certain  factors  that  encompass  the 
field  of  fracture  mechanics  and  their  use  in  achieving  the  capabili- 
ties of  Structural  Integrity. 

In  a flawed  or  cracked  material  an  evaluation  of  the  flaw  toler- 
ance (i.e.  fracture  toughness)  is  required.  An  analysis  of  stress 
state  for  the  loading  condition,  nature  of  crack  growth,  and  fracture 
criterion  must  be  considered  if  performance  or  life  prediction  are 
desired.  Stress  analysis,  and  particularly  stress  intensity  factors, 
are  evaluated  generally  by  either  analytical  and/or  numerical  methods 
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(Paris  et  al , 1965). 

After  obtaining  a suitable  representation  of  Kj , its  critical 
value  Kjc  must  be  determined.  Unlike  Griffith's  work,  where  a*/a 
and  /Ty  were  measured  independently,  Kjc  is  identified  for  most  engi- 
neering materials  by  a prescribed  fracture  test.  The  reason  is  that 
Gjc  cannot  be  measured  independently  as  Griffith  had  done  for  >. 

Obtaining  a material's  fracture  toughness  value  permits  a 
distinct  application  in  design.  If  design  stresses  can  be  accurately 
determined  for  the  structural  material  application  then  calculation 
of  the  largest  tolerable  flaw  size  can  be  established.  Inspection 
to  verify  that  no  flaws  larger  than  this  size  would  exist  frequently 
becomes  very  necessary  and  difficult.  Conversely,  the  minimum  in- 
spectable  flaw  size  can  be  used  to  set  the  allowable  load  limits 
for  the  structure.  It  would  be  ideal  if  flaws  in  structures  could 
be  detected  prior  to  ultimate  failure,  however,  this  is  not  always 
the  case.  An  arbitrary  division  of  the  three  toughness  regions  with 
their  associated  functional  failure  modes  are  cited  below: 


Failure  Mode 

Toughness  Region 

a. 

Overload,  plastic  yielding 

Plasti c 

b. 

Yielding  before  break 

Elasti c-Plasti c 

c. 

brittle  fracture 

Plane  Strain 

The  respective  consequences,  dependent  on  application,  could  bo: 
•i.  H.i inUMiur.oe , repair  and  overhaul 

b.  I n v i ronmental  contamination,  possible  structural  loss 

c.  Human  disaster,  catastrophic  failure 
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In  many  designs  the  probable  failure  mode  is  secondary  to  the 
number  of  cycles  to  failure,  however,  it  does  control  the  severity  cf 
failure  and  remains  the  primary  factor  in  assuring  failure-safe  de- 
sign. 

Crack  or  fracture  toughness  is  definitely  a material  property 
used  to  evaluate  structural  integrity  of  a material  providing  mode 
of  crack  growth  and  Griffith's  concept  is  applicable.  Due  to  incon- 
sistencies in  internal  microstructure , deviations  from  the  expected 
have  occurred.  Because  engineering  materials  are  often  non-glass 
like  and  Gjc  cannot  be  measured  independently,  much  work  has  been 
directed  to  repeating  Griffith's  analysis  for  materials  that  are  not 
perfectly  elastic,  brittle  bodies.  It  is  preferable  to  describe  frac- 
ture or  crack  toughness  in  terms  of  diverse  energy  spreading  mech- 
anisms with  accompanying  geometry  and  loading.  Progress  to  this  goal 
has  been  slow  due  to  limited  analytical  methods,  however,  many 
achievements  have  been  attained  due  to  precise  interpretation  of  ex- 
perimental results. 

Knowing  the  final  fracture  condition  with  its  attendant  para- 
meters is  important  in  identifying  the  mortality  of  a structure. 
However,  pursuing  a positive  aspect,  establishing  the  nature  of 
crack  growth  provides  information  with  respect  to  maintenance  cycles 
thus,  extending  time  to  component  wear  out  (Freudenthal .December 
1973  and  Hoeppner,  1973-1976).  This  study  on  17-4  PH  and  16-6  PH 
is  concerned  primarily  with  this  point  as  it  explores  crack  growth 
under  a loading  condition  termed  fatigue. 
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3.3  FATIGUE 

The  fatigue  failure  process  is  discussed  throughout  this  section 
to  develop  the  concept  of  fatigue-crack  growth.  Fatigue  failure 
occurs  when  a material  is  subjected  to  repeated  applied  cyclic  force. 
However,  fracturing  is  produced  by  a force-amplitude  that  is  far  be- 
low the  force  associated  with  fracture  under  a single  load  applica- 
tion. This  amplitude  decreases  with  increasing  number  of  repetitions. 
If  the  number  is  relatively  large,  fracture  occurs  well  within  the 
range  of  stresses  and  deformations  that  can  be  considered  elastic 
on  a "macroscopic"  scale.  These  fatigue  fractures  usually  occur 
without  noticeable  macroscopic  permanent  deformation  and  have  been 
produced  under  the  following  test  conditions: 

-fatigue  at  constant  load  amplitude, 

-fatigue  at  constant  strain  amplitude, 

-fatigue  under  spectrum  loading, 

-fatigue  at  positive  or  negative  mean  loads. 

Conventionally,  these  tests  have  been  conducted  with  notched  or  un- 
notched specimens  under  sinusoidally  constant  amplitude  load  cycles 
until  failure  occurred.  A curve  of  stress  versus  number  of  cycles  to 
failure  was  plotted  for  the  test  population  results. 

Recognizing  that  fatigue  from  design  discontinuities  and/or  in- 
herent material  or  nanufacturing  inconsistencies  progressively  damages 
the  material,  study  of  the  fatigue-crack  mode  was  undertaken.  The 
fatigue  crack  grows  incrementally  after  starting  from  a stress 
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concentration  and  spreads  gradually  by  mechanisms  of  slip,  cleavage, 
twinning  and  delamination.  When  the  remaining  cross  section  can  no 
longer  withstand  even  a single  further  application  of  the  cyclic 
force,  fracture  by  momentary  overload  begins.  It  is  at  this  stage 
that  the  crack  propagation  mechanism  associated  with  unstable  frac- 
ture under  a single  rapid  load  application  becomes  relevant  with  re- 
spect to  final  fatigue  fracture  (Freudenthal , 1958).  Here,  the  in- 
cremental crack  extension  ya  is  plotted  against  the  number  of  cycles 
(N)  to  yield  the  a-N  curve  (Figure  2)  which  then  can  be  differen- 
tiated to  establish  a crack  growth  rate:  Da/DU. 

Reflecting  on  section  3.2.2,  it  can  be  observed  how  the  study  of 
fatigue-crack  growth  relates  the  stress  intensity  parameter,  Kj,  to 
fatigue-crack  propagation  rates  until  fracture  where,  if  valid,  Kjc 
is  achieved. 

3.3.1  Fatigue-Crack  Propagation 

Crack  growth  processes  in  fatigue  can  be  expressed  as  a function 

of  the  fracture  mechanics  crack  tip  stress  intensity  factor,  K or  nK, 

versus  the  empirically  determined  crack  growth  rates  (Da/DN).  Com- 

Dlete  description  of  the  fatigue  stress  cycle  for  simple  wave  forms 

is  characterized  by  a maximum  stress,  S , and  a minimum  stress, 

Smin,  corresponding  to  a stress  range  AS  = (S[nax  - S|nin).  Now, 

associating  this  with  (K  „ - K.  ) = AK  permits  evaluation  of  the  ef- 
3 ' max  non  r 

fective  stress  intensity.  Additionally,  a stress  range  load  ratio 

S i 

(R  = ~ — ) can  be  computed  for  the  applied  stresses  and  used  for 
^max 


Figure  2 
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Crack  Growth  (Flaw  Size)  Versus  Number  of  Cycles  Illustrates 
Effect  of  Environment  and  Frequency.  Note  that  when  environ 
ment  is  influential  the  final  crack  size,  a , may  be  reduced 
(Hoeppner,  1975). 
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identification  of  crack  growth  rates.  For  zero-to-tension  cycling, 
R=0,  in  a non-aggressive  environment,  the  crack  propagation  generally 
follows  a typical  logarithmic  plotted  sigmoidal  curve  (Figure  3). 
Traditionally,  one  of  many  ways  the  second  region  (transition)  of 
the  curve  has  been  described  is  assuming  a linear  semi -empirical 
relationship  referred  to  as  the  power-law  or  relation. 

Da/ DM  = C(AK)m  (7) 

The  curve,  bounded  by  region  1,  illustrates  nonpropagating  fatigue 
cracks  when  Da/DN-*0  at  some  low  AK, . Above  region  3,  crack  growth 
accelerates  unstably  to  eventual  failure  at  some  high  AK^. 

The  parameters  C,  m,  AK.|  and  AK^  are  sensitive  to  the  material. 
The  exponent  (m)  generally  varies  approximately  between  2 and  5 for 
steels,  and  in  some  materials  is  as  high  as  lb  for  different  ma- 
terial conditions,  load  spectra  and  test  environments.  Generally, 
a small  value  of  (m)  leads  to  optimum  fatigue-crack  growth  resis- 
tance, although  there  are  definite  exceptions. 

The  value  (C)  is  considered  to  be  primarily  dependent  on  Young's 
Modulus  (E)  (Crooker,  1972).  However,  values  of  (C)  nay  vary  sub- 
stantially with  changes  in  alloy  groups  and  the  value  of  (m).  Often 
when  comparing  (Da/DN)  versus  (AK)  curves  of  two  Materials,  the  curves 
will  intersect  and  the  benefit  of  each  material  must  then  be  judged 
for  its  intended  application.  Consideration  should  be  given  to  the 
minimum  AK  values  that  apply  to  the  structural  situation  under  in- 
vestigation. Then  only  that  part  of  the  curve  above  the  minimum  :,K 


figure  3 Sigmoidal  Crack  Growth  Rate  Curve  for  Engineering  Materials 
or  Engineering  Alloys  Plotted  on  Log  Paper  (Crooker,  1973 
and  Kondas,  1974). 
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should  be  used  for  further  comparison.  This  type  of  approach  is 
important  in  highly  stressed  structures  where  the  flaw  detection 
capability  is  limited  and  resulting  AK  levels  are  in  a region  of 
medium  to  fast  crack  growth  rates. 

There  are  other  parameters  that  can  influence  the  shape  of  the 
curve  in  Figure  3.  Some  of  these  are  environment  (Figure  4),  stress 
ratio  (Figure  5),  and  load  wave  form  and  frequency.  In  general,  en- 
vironment has  its  greatest  effect  at  low  stress  intensities  and  low 
operational  frequencies  (Hoeppner,  1964,  1971  and  Crooker,  1972),  but 
exceptions  are  available.  There  are  many  contoi nations  of  the  above 
factors  with  effects  that  are  variable.  Some  of  the  combinations  will 
be  investigated  in  this  work. 

Presently,  to  determine  an  optimum  alloy  for  a fatigue  applica- 
tion, a materials  trade-off  may  be  conducted  using  conventional  para- 
meters such  as  n' , a Kjc,  oys  (Feltner  et  al,  1967). 
However,  this  approach  is  restrictive  because  it  lacks  the  necessary 
technique  to  optimize  or  indicate  those  parameters  that  could  be 
Improved  to  increase  the  fatigue  resistance  of  a new  metal 
(Grosskreutz,  1972).  This  must  be  done  by  examining.  Identifying  and 
relating  microstructures  to  specific  fatigue  resistance.  Frequently, 
when  we  have  attempted  to  improve  the  fatigue  resistance  by  raising 
the  fatigue  limit,  we  have  inadvertently  decreased  the  flaw  tolerance 
of  the  material.  In  addition,  the  material's  resistance  to  environ- 
mentally enhanced  subcritical  flaw  growth  has  decreased,  creating  a 
maintenance  problem. 
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Figure  4 Typical  Effect  of  Corrosive  Environment.  At  low  frequencies 
the  dashed  curve  shifts  to  the  left  and  generally  increases. 
Maximum  effect  is  at  low  AK  values  (Hoeppner,  1964-1971, 
Crooker,  1973  and  Kondas,  1974). 
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It  has  been  assumed  (Grosskreutz,  1972  and  Laird,  1967)  that  due 

to  high  stress  concentrations  present  under  high  crack  growth  rates 

-5  -5 

(i.e.  greater  than  4.0  x 10  in/cyclc  or  10.2  x 10  nm/cycle)  fa- 
tigue resistance  is  unaffected  by  microstructural  conditions  (i.e. 
second  phase  precipitates,  residual  stresses,  etc.). 

However,  to  alter  growth  rates,  inclusions  in  the  matrix  have 
been  observed  to  provide  initiation  sites  for  cracks  in  the  matrix 
(Harrigan,  1973  and  Hoeppner,  1969,  1975).  Thus,  if  the  inclusion  or 
interstitial  content  was  decreased,  leaving  a cleaner  alloy  matrix, 
the  fatigue-crack  growth  rates  should  decrease.  Further,  Hoeppner 
(1963)  has  maintained  that  the  grain  boundary  r.iisorientation  angles 
are  crack  direction  mechanisms  of  primary  influence  in  initiating 
flaws  by  controlling  the  energy  potential  to  its  least  value  thus, 
developing  slip,  cleavage  or  twinning  as  a result.  Outgrowth  of 
this  thinking  is  observed  in  areas  such  as  microstructural  and  mechan- 
ical texture  and  subgrain  boundary  cell  development  and  effect.  As 
expected  grain  size  has  been  observed  to  be  influential.  Grain  size 
effect  can  be  allied  with  slip  or  glide  mode  and  stacking  fault 
characteristics  as  presented  by  Hoeppner  (1963  and  1967)  and  Laird 
(1967).  If  the  metal  has  a high  stacking  fault  energy  (i.e.  wavy 
slip)  the  flow  stress  (static  or  dynamic)  is  essentially  independent 
of  grain  size.  It  appears  that  such  materials  form  a characteristic 
dislocation  cell  structure  that  depends  on  strain  amplitude  and  test 
temperature.  These  cells  are  much  smaller  than  the  grains.  Materials 
of  low  stacking  fault  energy  (i.e.  planar  slip)  show  a distinct  de- 
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pendence  of  fatigue-crack  growth  on  grain  size.  The  reason  is  that  in 
these  materials  dislocation  cells  do  not  form  at  high  cyclic  strains 
but  a uniformly  banded  structure  develops.  The  grain  size  effect  is 
considered  as  most  influential  for  planar  slip  mode  materials  such  as 
the  steels  being  examined  in  this  study  (i.e.  17-4  and  15-5).  How- 
ever, when  decreasing  grain  size,  fatigue-crack  growth  rates  decrease 
only  to  a point  depending  on  the  thermal  process  and  microstructural 
configuration  (Feltner  et  al,  1970). 

t 

! 3.3.2  Materials  and  Mechanics;  Orientation  and  Direction 

i 

To  achieve  an  increase  in  fatigue  life  by  decreasing  fatigue- 
crack  growth  rates,  an  optimum  approach  would  be  to  design  requisite 
metal  quality  through  microstructural  modifications  in  the  surface, 
grain  size,  phase  distribution,  texture  and  metal  cleanliness. 

Figure  6 graphically  diagrams  factors  to  be  reviewed  in  the  metal 

damage  process.  The  ability  to  relate  specific  microstructures  with 

I 

i fatigue  resistance  has  been  explored  minimally,  and  the  present 

i J 

\ 

' state-of-the-art  does  not  permit  accurate  predictions  of  the  above 

mentioned  for  optimum  fatigue  response.  This  field  will  have  to  be 
explored  by  asking  the  right  questions  and  utilizing  chemistry  and 

j physics  to  relate  to  exact  phenomenological  events  referred  to  for 
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Figure  6 Impacting  Factors  on  the  Physical  Process  of  Crack  Growth 
Under  Fatigue.  The  metal  quality  or  microstructure  is 
the  focal  point  from  where  fatigue  performance  diverges. 
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Chapter  4 


PROGRAM  PLAN 

4.1  TEST  PROGRAM  PURPOSE  AND  APPROACH 

The  design  ano  construction  of  Naval  Advanced  High  Performance 
ships  (AHPS)  presents  technological  similarities  to  the  airframe 
industry  where  minimum  weight  constraints  and  high  performance  re- 
quirements dictate  use  of  higher  strength  materials. 

Traditionally,  design  of  ship  structural  components  is  conser- 
vative to  keep  service  stresses  below  allowable  macroscopic 
material  property  limits,  such  as  ultimate  strength  and  0.2  percent 
yield  stress.  The  alloy  selection  process  sometimes  considers  other 
macroscopic  properties  such  as  fracture  toughness  and  ductility. 
Rarely,  however,  are  these  included  in  design  principally  because 
"traditional"  design  methodology  does  not  consider  these  quantities. 
The  "traditional"  engineering  design  practice  while  reasonable,  could 
be  improved  by  developing  a methodology  for  further  mechanical  prop- 
erty considerations. 

This  study  investigated  the  existence  of  a specific  parameter, 
dynamic  hold-time,  in  structural  material  applications  and  attempted 
to  identify  its  effect  and  micro-mechanism.  Generally,  the  hold- 
time parameter  is  not  considered,  but  it  has  been  found  to  be  a key 
factor  in  fai igue-crack  propagation  behavior  (Hoeppner  et  al , 1968, 
I'opp  i*l.  al,  l'H>9,  and  Ryder  <*t  al,  1973).  Unfortunately,  neither  an 
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individual  nor  combined  macroscopic  material  property  characteriza- 
tion of  conventional  design  properties  has  produced  evidence  to 
Indicate  hold-time  to  be  a culprit  In  numerous  service  condition 
cracking  Incidents  in  the  aviation  Industry.  However,  relatively 
simple  laboratory  scale  fatigue-crack  propagation  tests  can  reveal 
Interactions  between  applied  loads,  external  environment,  and  mic- 
rostructural  phenomena  and  help  forewarn  of  deterioration  In  crack 
growth  resistance  under  service  conditions.  The  generalized  concept 
for  use  of  fracture  mechanics  is  to  provide  a couple  between  idealized 
laboratory  and  in-service  fatigue-crack  propagation  behavior  to  enhance 
the  accuracy  of  life  prediction  procedures. 

The  objective  of  this  study  In  part  was  to  qualify  whether  existing 
circumstances  as  applied  to  candidate  alloys  of  hydrofoil  and  strut 
structures  produce  an  Increase  in  the  fatigue-crack  propagation  rate 
when  hold-times  typical  of  service  operation  are  considered. 

During  service  operation,  the  structures  experience  variations 
in  steady  loading  conditions.  In  the  hull  borne  mode,  the  foils 
and  struts  may  be  raised  or  lowered  for  a prolonged  time,  contingent 
on  ship  operation.  Thus  creating  an  increased  or  decreased  steady 
loading  on  a certain  area  of  the  structure.  Continuous  steady  loads 
can  develop  also  in  foil-borne  mode  through  biaxial  loading  and  drag 
bending.  If  cavitation  occurs,  a low  noise  ringing,  similar  to  hollow- 
blade  turbine  problems,  will  result  In  an  internal  friction  condition 
if  nit  properly  damped  (Bowie,  1974-1976). 
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To  simulate  this  type  of  load  behavior,  fatigue-crack  propaga- 
tion tests  were  conducted  with  various  hold-times  at  peak  loads. 

These  loads  and  stress  intensities  can  be  acquired  from  the  struc- 
tural analysis  of  the  hydrofoils  and  struts  providing  applicability 

i 

of  results.  In  stress  analysis  of  foil  and  strut  structures  con- 
tinuous perturbation  with  changes  In  temperature  gradients,  typical 
of  turbines,  is  not  a predominant  concern.  However,  in  view  of  fit- 
up,  welding  constraint,  heat  treatment  and  flame  straightening,  the 
stress  field  Is  not  necessarily  reduced  to  trivially  low  values  under 
a constant  temperature  and  at-rest  conditions.  Some  portions  of  the 
j structure  can  have  locked-ln  or  residual  stress  fields  of  magnitude 

\ approaching  tne  material  yield  strength.  Accordingly,  fatigue-crack 

l 

propagation  tests  will  be  performed  with  consideration  of  residual 
stress  variability  (R-ratlo)  as  a Impacting  feature.  The  R-ratlo 
will  be  patterned  after  that  in  the  relevant  and  applied  Navy  spectrum. 

4.2  DETAILED  DESCRIPTION 

* 

I 

4.2.1  Rationale 

I 

*. 

■ The  fatigue-crack  propagation  tests  were  conducted  In  two  modes. 

? First,  sine  wave  tests  performed  at  10  Hz  determined  base-line,  then  ! 

trapezoidal  wave  "hold-time"  test  were  conducted.  The  wave  form  dif- 
ference is  seen  in  Figures  7a  and  7b. 


Figure  7a  Sinusoidal  Wave  Form  Illustrating  P and  P . . 
* 3 max  mi  n 


Figure  7b  Trapezoidal  Wave  Form  with  One  Second  Rise  Time  and  One 
Minute  Mold-Time* 
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The  hold-time  and  stress  ratio  (R)  * e evaluated  In  accordance 
with  the  test  schedule  (Table  I).  This  s.udy  has  been  Integrated  with 
programs  at  Naval  Research  Laboratory  and  Lockheed  Rye  Canyon  Ad* 
vanced  Design  and  Research  Laboratory  to  provide  the  utmost  benefit 
to  the  U.S.  Navy  on  a Fleet  Support  level. 

The  dynamic  loading  procedure  used  In  testing  hold-time  enhance- 
ment of  fatigue-crack  propagation  should  not.  be  confused  with  static 
loading  procedures  of  sustained  load  (cracking)  flaw  growth  that, 
among  other  things,  consider  material  susceptibility  to  stress  cor- 
rosion cracking  (SCC)  at  a given  stress  Intensity  (Kjscc)*  Mechanisti- 
cally, the  two  events  may  be  similar  or  diametrically  opposite.  It 
has  been  demonstrated  that  alloy  response  to  stress  corrosion  crack- 
ing (SCC)  Is  sensitive  to  microstructure  (Krupp  and  Hoeppner,  1973). 
Whether  a microstructure  can  be  conditioned  to  prohibit  (SCC)  and 
mitigate  or  eliminate  hold-time  effect  In  fatigue  Is  an  area  of  on- 
going research. 

The  hold-time  effect  would  not  be  overly  Important  If  the  "Initia- 
tion of  cracks"  was  difficult  and  thus  a major  part  of  component  life 
spent  at  low  growth  rates  and  low  stress  Intensity  levels.  However, 

If  cracks  are  part  of  the  fabricated  structure  In  the  form  of  weld 
suck-back,  heat  treatment  cracks,  etc.,  then,  the  time  spent  In  "crack 
Initiation"  Is  lost  and  larger  initial  cracks  would  be  subject  to  re- 
latively high  stress  Intensity  levels.  Thus,  the  effect  of  hold-time 
would  become  quite  Important. 


TABLE  I TEST  MATRIX  FOR  EACH  MATERIAL 
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Multiple  cracking  frequently  observed  ahead  of  the  main  crack 
during  this  type  of  fatigue-crack  propagation  and  In  stress  corrosion 
tests  Indicates  an  Interesting  cracking  phenomena.  Such  multiple 
cracking  is  possibly  indicative  of  severe  multiple  cracking  under 
service  conditions  where  small  cracks  couple;  yielding  larger  cracks, 
higher  propagation  rates,  and  reduced  time  between  maintenance  periods. 

4.2.2  Materials  and  Heat  Treatments 

Metals  that  deform  plastically  by  slip  can  have  resistance  to 

j 

slippage  increased  by  any  treatment  that  will  harden  and  strengthen 
, the  metal.  There  are  numerous  ways  of  "hardening"  metals  such  as 

cold-working,  diffusion  hardening,  and  precipitation  hardening  (PH), 
j S^nce  PH  steels  are  the  subject  of  this  work  a discussion  of  the 

’ hardening  process  and  heat  treatment  process  used  on  them  follows. 

Precipitation  hardening  Is  caused  by  myriad  submicroscopic 
particles  being  deposited  on  slip  planes  throughout  the  alloy.  These 
* particles  strengthen  the  alloy  by  increasing  slip  resistance. 

I 

i Not  all  alloys  are  hardenable  by  precipitation  hardening.  Alloys 

1 

■ that  do  harden  by  this  mechanism  have  only  one  thing  In  common,  this 

Is,  a decreasing  solubility  for  one  phase  on  cooling. 

The  17-4  PH  and  15-5  PH  alloys  are  martensitic  PH  stainless  steels 
owing  their  properties  (Table  II)  to  a combination  of  two  Interacting 
metallurgical  phenomenon--the  transformation  of  austenite  (Pace-centered 

i 

Cubic)  to  martensite  (Body-centered  Tetragonal)  and  precipitation 

I 

J 
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hardening.  A precipitation  hardening  heat  treatment,  that  hardens 
by  precipitation  of  a copper-rich  phase  within  the  matrix  and 
tempers  the  matrix  structure  gives  these  alloys  their  properties. 

The  complete  heat-treating  sequence  consists  of  two  steps, 

viz; 

1.  Solution  Treatment  - 1900  F + 25  F (1038  C + 13.9  C) 
for  30  minutes,  then  air  cooTed. 

2.  Precipitation  Harden  - in  a range  from  900  F to  1150  F 
(482  C to  621  C)  for  1 to  4 hours,  then  air  cooled. 

The  microstructure  of  these  alloys  Is  controlled  by  chemical 
composition  and  heat  treatment.  The  constituents  that  comprise  the 
alloys  (Table  III  and  Figure  8)  can  be  classified  as  ferrite  formers, 
austenite  formers,  and  austenite  stabilizers  within  the  ranges  indi- 
cated In  Table  III.  The  balance  between  these  elements  Insures  a pre- 
dominant martensitic  microstructure  with  some  delta  ferrite  present  In 
17-4  PH  and  virtually  no  delta  ferrite  In  15-5  PH.  Carbon  and  nitrogen 
are  carefully  controlled  due  to  their  strong  Influence  on  the  auste- 
nite's stability  to  Insure  a full  martensitic  transformation  and 
develop  properties  to  their  fullest. 

At  the  solutlonlzing  temperature  (1900  + 25  F or  1038  C + 13.9  C) 
the  metallurgical  structure  of  17-4  PH  and  15-5  PH  Is  austenite  with 
dispersed  Islands  of  delta  ferrite.  The  solution  heat  treatment 
dissolves  copper,  the  hardening  element  In  the  austenite.  Upon  cool- 
ing from  this  temperature,  transformation  from  austenite  to  marten- 
site begins  at  approximately  270  F (133  C)  and  usually  Is  complete  by 
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TABLE  III.  MATERIAL  COMPOSITION 


ELEJEHIS 

17-4  PH* 

15-5  PH* 

VJFIGMT  PERCENT 

C 

0.07  MAX 

0.07  MAX 

fto 

1.00  MAX 

1.00  MAX 

P 

0.04  MAX 

0.04  MAX 

S 

0.03  MAX 

0.03  MAX 

Si 

1.00  MAX 

1.00  MAX 

Cr 

15.5  - 17.5 

14.0  - 15.5 

Mi 

3.0  - 5.0 

3.5  - 5.5 

Cu 

o 

in 

i 

o 

in 

2.5  -4.5 

Cb  + Ta 

0.15  - 0.45 

0.15  - 0.45 

Fe 

BALANCE 

BALANCE 

•vacuum-arc  remelted 


\ 

J 

I 

1 


I t 


lr  i if  1 litii  lit  k t hi.1  vN  ’ ij  -i. li«J-  i .Jdii 


WEIGHT  PERCENT 
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90  F (32.2  C).  Treatments  at  temperatures  below  1875  F (1024  C)  will 
result  in  reduced  tensile  and  yield  strength  properties  because  the 
copper  hardening  element  is  not  completely  dissolved.  Treatments  at 
tenperatures  above  1925  F (1052  C)  cause  excessive  grain  growth  re- 
sulting in  lower  tensile  ductility  and  reduced  impact  strength. 

After  solutionizing,  the  martensite  that  is  formed  on  cooling  is 
super-saturated  with  copper.  Subsequent  heating  to  temperatures 
below  the  precipitates  a copper-rich  phase  within  the  martensite 
increasing  hardness  and  strength.  The  stress  relieving  effect  due 
to  tempering  of  martensite  combined  with  precipitation  hardening  is 
important  in  terms  of  increasing  ductility  and  toughness  of  the  alloy. 
Heating  at  temperatures  above  causes  increasingly  more  au' 
to  form.  The  austenite  formed  at  these  temperatures  (1150-15-  t 
621  - 816  C)  does  not  transform  to  martensite  upon  cooling  due  to 
the  excessive  carbon  in  solution  stabilizing  the  austenite  (Armco 
Steel  Co.,  1969). 

4.2.3  Test  Specimen  and  Preparation 

Figure  9 shows  the  modified  compact  tension  specimen  used  in  this 
study.  It  is  referred  to  as  a wedge  open  loading  (WOL)  specimen. 

Note  that  this  specimen  is  slightly  longer  and  higher  than  the  com- 
pact tension  specimen  normally  used  for  Kjc  testing.  The  WOL  height 
to  width  ratio  is  H/W  = 0.486  versus  H/W  = 0.600  for  the  standard  Kjc 
specimen  (ASTM  E399-1972).  The  longer  specimen  (WOL  - ASTM  1975)  is 


Figure  9 Dimensioned  ledge  Open  Loading  (WOL)  Specimen  Used  in  Fa- 
tigue-Crack Growth  Testing. 
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preferred  for  crack  growth  testing  because  the  rate  of  change  in  Kj 
with  crack  length  is  not  as  rapid  as  in  smaller  specimens.  This 
slower  rate  of  change  makes  more  accurate  crack  growth  rate  measure- 
ments possible. 

All  test  specimen  blanks  were  cut  from  one  inch  (2.54  cm)  thick 
plates  such  that  the  direction  of  crack  propagation  is  parallel  to 
the  rolling  direction  (T-L  orientation).  In  addition,  a random  dis- 
tribution of  specimen  blanks  were  used  to  eliminate  possible  ordering 
effects  associated  with  specimen  location.  Figure  10  shows  the 
orientation  of  the  specimen  blanks  in  the  as-received  plate. 

After  solutionizing  all  test  specimen  blanks  within  an  argon 
atmosphere  the  blanks  were  final  machined  before  precipitation. harden- 
ing, again  under  argon.  A careful  examination  was  made  to  certify 
dimensional  control  on  these  speclnens. 

4.3  EXPERIMENTAL  TEST  EQUIPMENT 

4.3.1  Fatigue  Testing  Equipment 

A 20  kip  (88,960  nt)  closed  loop,  servo-controlled  electro  hy- 
draulic test  machine  was  used  for  these  tests.  The  machine  operated 
in  a load  control  mode.  An  amplitude  measurements  system  and  an 
oscilloscope  were  available  for  continuous  monitoring  of  load  during 
testing.  When  crack  opening  displacement  (COO)  was  monitored  the  load 
versus  COD  was  autographical ly  recorded.  The  output  signals  to  the 
recording  device  were  + 10  volts  full  scale,  single  amplitude. 


Figure  10  Fracture  Plane  Identification  where  L = Dii'ection  of  Grain 
Flow,  S = Short  Transverse  Grain  Direction,  and  T = Trans- 
verse Grain  Direction. 
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The  closed  loop  MTS  system  diagram  (Figure  11)  shows  the  load 
frame,  hydraulic  power  supply,  load  cell,  control  console  and  func- 
tion generator.  Note  that  the  oscilloscope  and  x-y  recorder  are 
included  in  the  test  assembly  (Figure  12). 

The  hydraulic  punping  unit  output  of  approximately  3 gpm  (11.4 
l/min)  under  pressure  of  3000  psi  (2072  nt/cm  ) is  used  to  drive  the 
hydraulic  actuator.  The  actuator,  similar  to  the  human  heart,  ac- 
cepts high  pressure  fluid  on  one  side  of  the  piston  actuator,  driv- 
ing it  in  the  opposite  direction.  The  control  of  the  actuator  is 
developed  through  a servovalve  and  a controller  that  compares  a command 
reference  signal  with  a feedback  signal.  After  both  are  satisfactorily 
summed  .nd  amplified  an  error  signal  is  generated,  inverted  and  sent 
to  the  servovalve,  metering  the  overall  flow  quantity  and  direction. 

The  entire  system  is  looped  closed  with  the  aid  of  feedback  signal 
from  the  load  cell.  Their  output  voltages  are  balanced  in  the  con- 
troller to  achieve  the  desired  response  profile. 
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4. 3. 1.1  load  Cell  Transducer 

The  force  applied  to  the  test  specimen  was  measured  with  a load 
cell  transducer  that  is  fatigue  rated  and  temperature  compensated. 

The  sensing  element  is  hermetically  sealed  within  the  load  cell  body 
and  consists  of  resistance  strain  gages  electrically  connected  to 
form  a balanced  Wheatstone  bridge.  The  signal  output  is  free  from 
the  effects  of  drift,  creep,  aspient  tenperature  changes,  hysteresis, 
or  resonance  from  statu  'uuJ'  •:  ip  to  SO  it;. 


i 
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Figure  11  Scaled  Drawing  of  Material  Test  System  Used  in  Compliance 
and  Fatigue- Crack  Growth  Testing. 
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4. 3. 1.2  Test  System  Calibration 

The  test  system  was  calibrated  prior  to  the  start  of  test  and 
at  four  month  intervals  during  testing.  Calibration  was  performed 
whenever  a transducer  or  test  system  component  Is  changed  or  anytime 
a testing  error  occurs  that  could  cause  the  system  to  exceed  trans- 
ducer limitations.  The  system  was  calibrated  as  a complete  unit 
from  the  point  of  input  to  the  transducer  to  final  signal  output  at 
the  recorder.  The  calibration  standards  that  are  used  are  traceable 
to  the  National  Bureau  of  Standards. 

The  load  and  the  COD  measuring  channels  were  calibrated  in  10 
equal  Increments  for  the  ranges  being  used.  Either  shunt  resistance 
calibration  techniques  or  master  load  cell  calibration  techniques  were 
used  to  calibrate  the  load  weighing  system.  COD  system  calibration 
was  performed  with  a frame  micrometer  with  0.0001  inch  (0.000254  cm) 
or  less  thimble  graduations.  The  maximum  allowable  error  at  the  test 
system  recorder  was  +0.5%  of  full  scale  single  amplitude  at  each  of 
the  10  increments.  The  COD  gage  was  calibrated  at  the  temperature 
and  the  environment  in  which  the  testing  was  performed.  The  load  cell 
was  adequately  protected  to  ensure  that  ambient  conditions  exist  at 
the  load  cell  during  all  testing. 

4.3.2  Accessories  for  Fatigue  Testing 
4.3.2. 1 Clevises 

The  load  was  applied  to  the  test  specimen  with  two  > atched  cle- 
vises. The  load  line  went  through  the  center  of  each  clevis.  The 
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loading  holes  were  0.50  inch  (1.27  cm)  diameter  with  load  pins  not 
less  than  0.484  inch  (1.23  cm)  diameter.  The  clevis  was  capable  of 
maintaining  specimen  vertical  alignment  witii  spacers  between  the  . 
specimens  and  clevis  inside  walls. 

4. 3.2.2  Crack  Opening  Displacement  Gage 

The  edge  crack  opening  displacement  (COD)  was  measured  with  a 
double-cantilever,  resistance  strain  gaged,  clip-on  displacement  gage 
transducer.  The  gage  was  similar  to  the  displacement  gage  described 
in  reference  (ASTM  £399-1972).  However,  during  ttie  data  reduction  the 
thickness  of  the  knife  edges  used  to  mount  the  COD  gage  was  taken  into 
account. 

4. 3.2. 3 Gaertner  Microscope 

This  microscope  was  used  at  33x  magni fi cation  and  mounted  in  a 
micrometer  slide  for  optical  recording  of  crack  extension.  Due  to  the 
motion  of  the  system  and  the  specimens,  scribe  marks  were  placed  at 
every  0.100  inch  (0.254  cm)  from  the  crack  tip  to  crack  run  length. 
These  markings  assist  the  crack  reading  accuracy  because  of  decreased 
transit  from  a known  point  to  the  crack  tip.  The  definition  of  crack 
tip  location  is  known  to  be  a function  of  crack  tip  displacement,  test 
material,  condition  of  specimen  surface  and  lighting  source.  For  this 
reason  these  factors  were  of  primary  concern  during  testing. 
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4.4  TEST  PROCEDURES 

All  test  specimens  were  precracked  under  sinusoidal  cyclic  load- 
ing prior  to  generating  the  crack  growth  data.  The  maximum  stress  in-' 
tensity  factor  during  precracking  did  not  exceed  17  ksi  /TrT  (18.7 
MPa  rm)  with  the  corresponding  stress  ratio  (R)  equal  to  0.05.  The 
test  frequency  used  for  precracking  was  10-15  Hz.  The  modified  com- 
pact tension  specimens  were  precracked  approximately  0.050  to  0.075 
inches  (0.127  to  0.190  cm)  beyond  the  end  of  the  machined  notch. 

After  precracking  to  an  acceptable  length,  a variation  in  the 
cyclic  loading  or  frequency  was  required  to  produce  a beach  mark  at 
the  precrack  front  to  evaluate  this  crack  front  eccentricity  or  bowing 
shape  after  final  fracture. 

following  precracking  and  marking  the  crack  front,  fatigue- 
crack  growth  testing  was  conducted.  Crack  growth  data  (a  vs.  N) 
was  recorded  optically  for  both  back  and  front  crack  extensions. 

Crack  length  measurements  were  made  at  crack  growth  intervals  not 
exceeding  0.01  inch  (0.025  cm).  Both  sides  of  crack  length  did  not 
have  a difference  exceeding  0.05  Inch  (0.127  cm)  at  any  time  during 
testing. 

For  testing  in  a salt  water  environment,  the  specimens  were 
placed  in  a sealed  plastic  bag  or  for  horizontal  load  frames  the 
specimens  were  dipped  into  the  salt  water  tank  (Figure  13). 

The  test  region  of  the  specimen  was  immersed  in  a solution  of 
3.5*  NaCl  in  deionized  water  which  was  maintained  at  a constant  flow- 


Figure  13  The  Portion  of  the  Specimen  to  be  Exposed  to  the  Salt  Water 
Environment  is  Illustrated  in  Position  within  the  Tank. 

The  water  is  fed  through  and  drained  without  reuse. 
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through  rate  of  approximately  one  half  pint  per  hour  (0.214  1/hr). 

The  solution  was  not  reused. 

Test  Interruptions  were  avoided  where  practical,  since  transient 
periods  of  crack  growth  retardation  have  been  observed  following  these 
Interruptions. 

The  maximum  and  minimum  load  excursions  were  monitored  con- 
tinuously during  the  test  by  a memory  circuit.  This  was  able  to 
detect  and  record  load  spikes  which  might  arise  «rom  extraneous 
electrical  or  mechanical  sources  even  though  limit  detectors  will 
monitor  gross  deviations  for  shut-down. 

During  the  test-  specimen  compliance  was  developed  by  measur- 
ing crack  opening  displacement,  COD.  The  process  will  be  discussed 
more  thoroughly  In  chapter  5. 

The  mathematical  relationships  for  the  specimen’s  stress  Intensity 
calibration  will  be  explored  In  chapter  5 as  twu  types  of  approaches 
are  considered.  The  Wilson  polynomial  and  Weibull  exponential  equa- 
tion* comprise  the  two  data  reduction  techniques  considered  for  this 
work.  The  rationale  for  these  techniques  will  be  discussed  in  the 
following  section. 
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EXPERIMENTAL  OATA  REDUCTION  DESIGN  AND  PROCEDURES 

5,1  WE  I BULL ' S SURVIVORSHIP  FUNCTION 

Deviations  of  the  Weibull  survivorship  function,  termed  the 
third  asymptotic  distribution  of  smallest  values,  are  given  in 
Gumbel's  text  (1958).  Weibull's  classic  paper  in  1949  was  the  pro- 
genitor for  present  applications  to  statistical  representation  of 
fatigue  failures.  Analytical  procedures  for  estimating  the  three 
parameters  of  the  Weibull  survivorship  function  were  sunwarized 
in  1963  (ASTM  1963). 

The  rationale  as  to  why  the  Weibull  function  has  been  pursued 
may  be  obvious  with  a succinct  overview  of  existing  expressions  and 
their  inherent  difficulties.  The  most  common  is  attributed  to  Paris 
(1965) 

Da/DN  = C(AK)m  (8) 

where  C and  m are  acquired  from  the  data  plot.  Here,  a significant 
deviation  occurs  as  the  data  does  not  conform  to  the  equation  at  low 
and  high  Da/L)N  values.  It  also  tends  to  mask  effects  due  tc  material, 
frequency,  envi ronment, 1 oad  spectrum,  and  stress  ratio.  The  most 
common  formulation  for  stress  ratio  correction  used  for  Da/UN  Inis 
been  proposed  by  Forman  (1967): 
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Da/Ofl  = k (9) 

c " max 

Even  though  applicable  at  higher  crack  growth  rates,  it  tends  to 
be  in  error  at  low  cracking  rates  where  the  largest  inpact  occurs 
to  long  life  structure  (Bnwie  et  al,  1974).  Another  expression 
for  stress  ratio  effect  was  explored  by  Walker  (1970): 

1 — in  n [n.  ni, 

40  ■ ■ ‘W’-R>  CO) 

This  approach  gives  a good  fit  to  the  data  at  low  crack  growth 

rates.  However,  as  K „ increases  to  a critical  value  and  environ- 

max 

mentally  assisted  cracking  effects  are  considered,  the  equations  will 
be  in  error  (Jaske  et  al , 1973). 

Overall,  a number  of  curve  fitting  numerical  relations  have  been 
proposed  to  assist  the  graphical  mechanics  of  fatigue-crack  propaga- 
tion. However,  no  single  relation  has  gained  acceptance  by  the  fa- 
tigue and  fracture  mechanics  conmunity. 

The  fitting  power  of  Weibull's  relation  is  demonstrated  by  ob- 
serving the  weaknesses  in  the  data  fitting  for  traditional  poly- 
nomial stress  intensity  calibrations.  This  evaluation,  will  be 
explored  in  the  following  section  using  ompliance  measurements. 
Additionally,  benefits  and  physical  rationale  are  suggested  for  the 
fatigue-crack  growth  life  predictions  approach. 
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5.2  COMPLIAHCE  MEASUREMEI/T 

5.2.1  Theoretical  Derivation  of  Conpliance 

An  entirely  new  relationship  for  the  representation  of  compli- 
» 

ance  measurements  is  considered  herein.  The  relationship  incorporates 
a Wei  bull  type  function  with  the  associated  Weibull  parameters  k,  e, 
and  v.  Before  discussing  the  approach,  a brief  review  of  the  histori- 
cal mathematics  for  compliance  calibration  to  supplement  section  3.2.1 
will  be  given.  Examining  specific  erergistic  relationships  of  the 
b-concept  (i.e.  crack  extension  energy),  compliance  determination  may 
be  considereo.  Assuming  a center  cracked  panel  or  a single-edge 
notch  configuration  (Figure  14a)  where: 

c = 2a  total  crack  length 
6c  crack  extension 
eQ  elastic  elongation 

Tne  increment  of  energy  consumed  in  crack  extension,  6c,  is  the 
change  of  strain  energy  with  crack  extension  for  fixed  grip  condi- 
tions. Plotting  load  (P)  versus  displacement  (eQ)  (Figure  14b)  re- 
sults in  an  energy  relation 

u 1/2  Peo  (11) 

Approaching  the  problem  by  utilizing  a spring  constant.  M,  beyond 
the  free  length  to  insure  linear  response  yields: 
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P = MeA 
0 

(12) 

P/M  = eQ  constant 

(13) 

The  compliance  of  a specimen  is  designated  as  1/M  in  this  approach. 

Differentiating  equations  (11)  and  (13)  with  respect  to  (c)  for 
a fixed  (e)  yields 


where  = 0 

Substituting  equation  (13)  into  (14) 


04) 


(IB) 


(16) 


(17) 


Recognizing  the  equivalence  of  terms  in  equation  (17)  to  equation 
(lb) 


(18) 
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The  parameters  P and  M are  defined  per  unit  thickness.  A plot  of 
specific  loads  and  displacements  is  shown  in  Figure  15. 

By  using  the  compliance,  (1/M),  and  plotting  it  against  the 
crack  size  (i.e.  c^ , i = 0,1, 2, 3),  the  failure  parameter,  Gc  for  a 
flawed  material  can  be  graphically  illustrated  in  Figure  16.  Here, 
G is  the  crack  extension  or  driving  energy,  where  G = - ~ ; u = 

p 

energy  internal  to  body. 


Thus,  the  critical  crack  driving  energy  can  be  represented  as 
f ol 1 ows , 


5.2.2  Experimental  Analysis  in  Weibull  Parameterics 

The  future  evaluation  in  crack  extension  energy,  stress  intensity, 
crack  growth  and  fatigue  life  predictions  has  a critical  focal  point 
at  compliance  de term-' nation.  Traditionally,  the  polynomial  fit  has 
been  employed  to  a moderate  degree  of  accuracy  over  a specific  range 
of  crack  growth.  An  attempt  to  increase  the  precision  of  the  curve 
fit  and  expand  the  range  led  to  the  following  work.  The  rationale  to 
pursue  a Weibull  type  function  stems  from  the  empirically  oriented 
concept  of  preference  for  best  fit  to  date  with  a reduced  or  minimum 
number  of  fitting  paraiieters.  A purely  mathematical  advantage  of 
the  Weibull  type  function  is  observed  as  all  data  curves  can  be 


Figure  15  Compliance  Illustration  of  Increasing  Crack  Length  (CQ  to 
C3)  Versus  Load-Displacement  Relationship. 


Figure  16  Compliance,  the  Reciprocal  of  Stiffness,  Fitted  to  a Tra- 
ditional Polynomial  over  Increasing  Crack  Size. 
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plotted  through  a coronon  point  permitting  future  relationships  to  be 
dependent  on  strictly  one  parameter,  shape. 

The  algebraic  mathematics  appear  rigorous  compared  to  a series 
expansion  of  the  polynomial  type  function,  but  is  relatively  elemen- 
tary in  the  manipulations  for  compliance  determination. 


Base  equation:  1 - a/W  = Exp 


loge{loge( 


EB6 


+ e)} 


(21) 


The  first  order  is  to  establish  through  linear  regression 
analysis  a coefficient  of  correlation,  r,  for  sample  of  data.  In 
this  case  the  sample,  N,  was  26  data  points  measured  by  J.T.  Ryder, 
Lockheed-Cal i fornia  Rye  Canyon  Advanced  Design  and  Research  Labora- 
tory. 

The  base  equation  (21)  was  modified  to  reflect  the  below  form 
(22)  with  variables  of  the  regression  Y,  yy  indicated 


loge(l 


a/W) 


1/k 


= 1/v 


’lo9e,lo9e(y 


(22) 


Y = YY 

Y = bX  + a 


The  computations  and  data  insert  are  carried  out  through  a re- 
mote IBM-2741  terminal.  The  dialogue  is  in  the  Time  Sharing  Opera- 
tion language,  a TSO  version  of  PL/I.  The  example  of  Weibull  para- 
meter estimation  employs  one  stored  subroutine. 
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As  the  program  is  called  and  entered,  the  number  of  sample 
points,  H,  are  introduced.  The  option  for  center  crack  panel  or 
wedge  open  loading  specimen  is  presented,  WOL  is  selected.  After 
corresponding  (c./W)  versus  (CEB)  (i.e.  EB6/P)  values  are  input  to 
the  system  the  iterative  loop  displays  numerous  values  of  w,  which 
is  1/k,  versus  r.  For  structural  alloys  studied  to  date,  the 
Weibull  shape  parameter  k = 1/w  shows  0-w-l. 

Upon  reaching  the  best  w and  r,  the  cycle  is  terminated  and  the 
best  w is  supplied.  At  this  point  the  regression  variable  forms  are 
printed  along  with  26  values  of  Y and  YY.  The  values  are  plotted  in 
Figure  17  and  show  only  one  point  outside  the  scatterband,  1.2  is 
permitted. 

The  following  pages  contain  the  slope  of  regression  line, 
scatter  in  the  vertical  direction  of  observed  points  about  the  re- 
gression line  linear  equation  and  its  measure  of  fit  and  WOL  para- 
meters k and  v. 

Through  all  remaining  plots  the  k and  v will  be  incorporated. 

The  first  demonstration  of  this  comes  when  the  pleasured  CEB,  a/W, 
and  k and  v terms  are  implemented  to  return  a calculated  a/W.  The 
calculated  values  were  incremented  over  an  entire  range  of  a/W, 
designated  normalized  crack  length,  and  plotted  in  figures  IB  and 
19  in  conjunction  with  22  of  the  data  points  from  the  measured  sample. 
The  correlation  is  impressively  precise.  Knowing  that  specimen  com- 
pliance is  the  integral  part  of  a stress  intensity  calculation, 
three  compliance  fitting  equations  are  equalized  and  compared  as 
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figure  17  Regression  Variables  Graphed  with  Correspcndino  26  Data 
Points  and  a 2<i  Scatter  band. 
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Figure  18  normalized  Plot  of  Compliance  Versus  Crack  Length.  The 
circled  11  points  characterize  11  measurements  performed 
by  J.T.  Ryder  (Personal  Comnunication) . 
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figure  19  Normalized  Plot  of  Compliance  Versus  Crack  Length  over  an 
Extended  Range.  The  circled  11  points  were  performed  by 
J.T.  Ryder  (Personal  Communication). 
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5.2.3  Summary  on  Stress  Intensity  Calibration  factor 


In  showing  the  entire  range  of  data  in  Figure  20  and  to  demon- 
strate the  mathematical  advantages,  the  regression  variable  YY  is 
plotted  against  the  normalized  crack  length,  a/W.  Because  at  a/W, 
1-1/e  = 0.6321  where  e,  here  is  the  naperian  base  equal  to  2.71 
while  YY  is  equal  to  1.0,  all  curves  must  pass  through  this  common 
point  leaving  k,  the  Weibull  shape  parameter  variable.  Addi- 
tionally, the  exponential  approach  shows  a higher  degree  of  cor- 
relation to  data  fit. 

Expanding  on  this  final  item,  it  was  determined  that  a compari- 
son should  be  developed  between  the  values  designated  for  two 
heavily  used  polynomials  and  the  exponential  expression  of  Weibull. 
The  two  expressions  are  *2  and  h 3 of  figures  21  and  22.  The  later 
will  be  the  polynomial  for  the  WOL  specimen  used  in  this  study.  The 
numerical  relations  of  the  curves  are: 


«1) 


a dCEB  |1/Z 
W d a/W 


CEB  = EB*/P 

= Exp  {Exp( v[-loge(l-a/W)]1/k) ) - Exp(l) 
where  v = 1.46158123 
k - 2.99580674 


"2 ) 


v3  W 

where  Y 29.6 
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Figure  20  Regression  Variable  Versus  Crack  Length  Plots  all  26  Points 
Over  the  Entire  Range  of  a/W. 


Figure  21  Comparison  of  Mathematical  Compliance  through  Numerical 
Relations  of  C3  Curves  Plotted  over  a Normalized  Crack 
Length  Ronge  0.1  to  0.6. 


Figure  22  Comparison  of  Mathematical  Compliance  through  Numerical 
Relations  of  C3  Curves  Plotted  over  a Normalized  Crack 
Length  Range  0.5  to  1.0. 
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*3)  C3  = S ' Y (Wilson  Polynomial) 

where  Y = 30.96  - 195.8  (~)  + 730.6  (g)2  - 1106.3  (jy)3 
+ 754.6  (|)4 


These  polynomials  (#2  and  #3)  are  empirical  best  fits  to  data  repre- 
senting different  specimen  dimensions. 

The  rationale  for  multiplication  of  the  polynomial  values  Y and 
a/W  is  as  follows  from  Irwin  and  Kies,  1954 


This  is  often  c , essed  as: 


v . ..  BW  .. 

i “ or  K. 

1 1 


Thus,  it  is  possible  to  write 


Y 


(a/W^ 


(24) 


(25) 


Thus,  if  oO  then  Y Knowing  this  and  wanting  to  erradicatc  singu- 
larity, was  used  for  comparison.  (Mote:  Analytic  functions  may 

have  different  types  of  singularities.  A singular  point  of  an  ana- 
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lytic  function  f(x)  is  a point  where  f(x)  ceases  to  be  analytic  such 
as  non-di fferentiable.  Also,  f(x)  is  singular  or  has  singularity  at 
that  point.  The  function  f{x)  is  termed  singular  at  infinity  if 
f(-jj-)  is  singular  at  x = 0).  Taking  advantage  of  Wessel's  work  with 
the  Wilson  polynomial  and  combining  with  the  modification  in  equa- 
tion (23)  an  expression  for  can  be  written  as 


C3  « (a/W)  -Y  (27) 


The  polynomial  equations  were  modified  as  stated  to  provide 
equivalence  for  comparison  and  plotted  in  Figures  21  and  22.  The 
stress  intensity  calibration  parameter  is  expressed  as 
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The  error  or  difference  plot  in  Figure  23  shows  the  most  exten- 
sive deviation  of  the  polynomials  outside  a/W  range  of  0.3  to  0.7 
where  polynomials  cannot  account  for  the  response.  A large  difference 
is  noted  between  curves  1 and  3 where  3 is  the  new  polynomial  ex- 
pression moving  toward  adoption. 

t 

5.3  LIFE  PREDICTION:  REFINING  NUMERICAL  ANALYSIS 

5.3.1  Convergence  on  Fitting  Quality 

5. 3. 1.1  Introduction 

This  section  provides  a review  of  some  numerical  interactions 
regarding  fitting  relations  and  their  quality.  By  standardizing  the 
discrimination  power  through  precise  data  reduction  and  evaluation 
techniques  and  increasing  measurement  efficiency,  it  is  considered 
possible  to  provide  a methodology  to  reliably  probe  the  macro-mcro 
interface  of  materials  properties  with  respect  to  their  fatigue-crack 
growth  behavior. 

5. 3. 1.2  Candidate  Exponential  Curve  Fitting  Relations 

It  is  usual  to  plot  fatigue-crack  propagation  test  results  in  the 

form 

Da/DN  » f (\K)  (32) 

where  Da/l)U  = crack  propagation  rate  and  :»K  = stress  intensity 
range.  A suitable  curve  fitting  relation  for  this  type  application 
is  one  with  a mathematical  fit  to  the  total  physical  behavior  of  fa- 


Figure  23  Percent  Differences  Between  the  empirical  Calibr^ofT Para- 
meters for  Stress  Intensity  over  a Normalized  Crack  Length 
Range. 
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tigue-crack  growth.  This  relation  must  have  the  following  properties: 

(1)  The  relation  should  lend  itself  to  distinction  be- 
tween slow  and  fast  crack  growth  rate  segments. 

(2)  As  AK  approaches  an  upper  limit,  1^,  the  crack  growth 
rate  should  increase  indefinitely. 


1 im 

^ > 1 

§/  ' 


Da 

DM 


(33) 


(3)  As  AK  approaches  a lower  threshold  limit,  Kth,  the 
crack  growth  rate  should  diminish  to  zero 


1 im 
AK  -K 


Da 

DM 


- 0 


th 


(34) 


A new  class  of  exponential  curve  fitting  relations  with  mathe- 
matical form  in  the  fatigue  data  reduction  Weibull  distribution 
function  were  explored  (Bowie,  1974-1976).  This  form  is  the  same  form 
similar  to  that  discussed  for  compliance  considerations.  Three  mem- 
bers of  the  class  are: 


Type  Relations 

A)  1-F  = Exp  , (36) 

or  II  = e + (v  - e)  |(  - loge(l  - f))  | 
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f / log  (H  +1)  - e \ k 
B)  1-F  ■ Exp  - 6 „ _ . ) 


v - e 


or  H = -1  + Exp  (e  + (v  - e)  |(-  1 oge( 1 - F))  ) (36) 


C)  1-F  = Exp 


1o9el°9e(H  + exp(l))  - e \k 

v - e 


or  H * -Exp(l)  + ExpExp  (e  + (v  - e)  (-  1 oge ( 1 - F))  ) 

These  relations  meet  the  three  criteria  of  acceptability  previously 
Stated.  Computations  are  performed  at  an  IBM  Type  2741  remote  com- 
puter terminal  in  PL/I  language.  The  fitting  routine  contains  logic 
for  estimating  the  curve  fitting  parameters  k,  e,  and  v by  a combina- 
tion of  linear  regression  and  correlation  coefficient  optimization. 

In  practice,  Da/DN  and  AK  values  of  a data  set  are  stored  as  one 
dimensional  arrays  DaDN  and  dK.  The  regression  analysis  is  per- 
formed on  variables  X and  Y,  in  terms  of  the  equation  of  regression 


( - loge  Y)w  = bX  + a , 


which  has  alternative  form 


Y = Exp  (-  (^-|)k)  (39) 

where  k = 1/w, 

e = -a/b, 

and  v = (1  + be)/b 

The  following  operations  are  performed,  for  example,  when  applying 


candidate  Type-B  curve  fitting: 
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DO  i > 1 TO  H 

X = loge(DaDN(i)  + 1) 

Y = (-  loge(l  - dK(i)/Kb))w 


Regression  and  Correlation  Coefficient  Logic 


END 


In  this  case,  the  routine  yields  the  Type-B  relation 


1 - AK/Kfa  = Exp 


' log  (Da/DN  + 1)  - e 


v - e 


(40) 


where  = F,  and  Da/DN  - H 


Parameter  is  the  stress  intensity  range  where  Da/DII  is  indefinitely 
large.  Estimation  of  requires  exercise  of  judgement  in  addition 
to  interactive  computing  effort.  At  the  outset  of  data  set  analysis, 
all  three  types  of  exponential  curve  fitting  relations  are  tested, 
with  F = AK/K^  and  H = Da/DN  for  each  type.  Measures  of  goodness  of 
fit  are  used  i r.  deciding  which,  if  any,  type  is  suitable.  For  all 
three  types,  the  particular  stress  intensity  range  where  Da/DN  = 0 
is  called  the  threshold  stress  intensity,  K^,  where 


vth 


= K.  (1  - Exp  (-  (- 


-e 


'))  ) 


v - e 


(41) 
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Physical  significance  is  attached  to  threshold  K^h  for  cases  where 
threshold  parameter  e _<  0. 


a. 

For  e 

= 0’  Kth 

= 0 

(42a) 

b. 

For  e 

< 0,  Kth 

> 0 

(42b) 

c. 

For  e 

> 0,  the 

suggestion  is  that 

(42c) 

crack  growth  begins  at  some  Da/DN  > 0 
for  arbitrarily  small  AK  > 0. 


5. 3. 1.3  Measures  of  Goodness  of  Fit 


For  a data  set  of  sample  size  Np,  linear  regression  is  performed 
with  systematic  iteration  on  parameter  w,  (w  = 1/k),  until  a sample 
correlation  coefficient,  r is  maximized.  In  the  routine 


N 


r = 


Np 

P J,  V, 


II 


P 

l X, 


i = l 


i = 1 


Y. 

l 


NP 

J,  V 


Nr 


HP  . 


Np 

( l 

i =1 


V 


(1/2 

W) 


To  the  extend  that  it  might  be  reasonable  to  assume  that  Y = 

1 /k 

[ - 1 oge ( 1 - AK/K^)]  ' for  optimum  k is  randomly  distributed,  it  may 

2 

be  inferred  that  regression  of  Y on  X accounts  for  100  r percent  of 
the  variance  of  Y.  The  scatter  of  Y . , i = 1,2  ...  Np,  about  re- 
gression line 

Y = bX  + a 


is  measured  by  computing  the  variance  term 
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2 1 
“2  * fa  <’ 


♦ v V11' 


(44) 


It  is  a good  idea  to  count  the  number  of  points  (X^,  which  fall 
between  the  lines 

Y*  = bX  + a + 2o 
and  Y"  = bX  + a - 2o 


When  approximately,  95%  of  the  points  do  fall  between  lines  Y*  and  Y", 
it  is  reasonable  to  plot  2a  bands  in  plots  of  Da/DN  versus  AK.  For 
example,  Type-B  curve  fitting  yields  2o  bands  defined  by  the  relations 


(Da/DN)'  = -1  + Exp  (e*  + (v 
and  (Da/DN)"  = -1  + Exp  (e"  + (v 

where  e'  = -(a  + 2o)/b 

v'  « (1  + b e')/b 
and  e"  = -(a  - 2o)/b 

v"  = (1  + b e")/b 


- e')  ( -loge(l  - AK/Kb))  1/k) 

- e")  ( -loge(l  - AK/Kb) ) 1/k)  , 

(45) 


the  aid  of  the 
(Yi  • (b  X.  + 


A third  measure  of  goodness  of  fit  is  monitored  with 

routine.  It  is  the  s 

/* 

a) ) , and  is  called  0£ 


routine.  It  is  the  sum  of  squares  of  the  deviations 
,2 


Np 

i.e.,  D2  = 2T  (Y.  - (b  X.  + a))2 
i-l  1 1 


(46) 


For  Type-B  curve  fitting 


Other  investigators  may  choose  to  adopt  an  alternative  computational 

procedure  designed  for  selection  of  best  values  of  parameters  Kb»  k, 

2 

e,  and  v by  minimizing  0 . 


5. 3. 1.4  Distinction  Between  Slow  and  Fast  Crack  Growth 

In  current  damage  tolerant  design  practice,  the  terms  "slow 
crack  growth"  and  "fast  crack  growth"  often  are  used  in  an  entirely 
descriptive  way.  The  following  quantitative  procedure  represents  a 
significant  analysis  development  for  distinguishing  between  slow  and 
fast  growth.  The  procedure  must  be  applied  with  caution,  because  re- 
sults are  highly  dependent  on  definition  of  the  parameter  the 
stress  intensity  range  where  Da/DN  increases  without  bound.  A given 
data  set  may  or  may  not  contain  sufficient  information  to  warrant  se- 
lection of  plausibly  physical  value  K^.  In  the  event  that  acceptable 
Kb  does  not  appear  to  be  intrinsic  to  a data  set,  it  is  necessary  to 
select  Kb  arbitrarily.  It  is  believed  that  possible  difficulty  in 
estimating  an  unbiased  value  has  a parallel  in  ASTM  C 399 , where 
it  is  necessary  at  times  to  specify  that  a fracture  mechanics  test  ha; 
yielded  an  invalid  estimate  of  fracture  toughness  KIc*  The  way  in 
which  a value  K.  is  selected  when  applying  the  routine  to  a data  set 
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( DaDN ^ , dK. ) is  as  follows:  The  largest  value  of  AK  in  the  set  is 

noted  and  an  initial  estimate  for  K.  is  chosen  to  be 

b 

K.  = AK  , + c,  0.01  < c < 1.0  (48) 

b max  — — 

An  initial  curve  fitting  relation  is  found,  and  the  initial  correla- 
tion coefficient,  r , is  noted.  The  computational  process  is  repeated 
with  successively  larger  estimates,  usually  in  increments  of  10 
ksi  /Tn  (11  MPa  /ni).  Two  possibilities  exist: 

* Either  r can  be  maximized  at  a particular 

• r increases  relatively  slowly  in  a monotonic  way,  even 

for  values  K,  >>  AK 

b max 

In  the  first  instance,  is  optimized  to  + 0.5  ksi  /in  (^0.55  MPa 
ri)  approximately.  In  the  second,  appeal  is  made  to  similar  data  set 
examinations,  to  the  approximation 

= H - Load  Ratio,  R)  K,  (49) 

where  K might  be  a fracture  toughness  value  (e.g.  Kjc  or  K^)  deter- 
mined in  some  other  test,  or  is  selected  arbitrarily,  sometimes 
with  a yield  stress  value  in  mind. 

b.3.1 .b  Risk  Analysis 

The  ratio  of  highest  to  lowest  Da/DN  value  in  a data  set  can  be 
of  the  order  of  10^.  It  is  no  wonder  that  fracture  mechanics'  in- 
vestigators plot  Da/Dfl  to  a logarithmic  scale.  Whether  a linear  or 
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logarithmic  ,'.K  scale  is  used  is  a matter  of  custom.  The  log-log 
format  can  be  helpful  when  the  data  seem  to  be  reduced  to  a straight 
line  o/er  a control  /-.K  region.  The  log-linear  format  is  recomnended 
when  exponential  curve  fitting  is  appropriate. 

The  trend  of  the  natural  logarithm  of  Da /Oil  witt  increasing  ,.K 
is  used  as  a guide  in  defining  the  transition  between  slow  and  fast 
growth.  The  trend  is  called  a risk  function,  r*,  given  by 


* 4 loye  Da/L)U  1 (IJa/DIl) 

' ' ' ' Tba/oiiy 

For  Type-B  curve  fitting. 


(bO) 


oj  (Jto/DU  * JJ 
k Ua/DN 


""0"- 


.-1  * 1/k 


'K'V 

Tk/'k^T" 


(bi) 


When  plotted  as  a function  of  /.K/K^  or  ..K,  risk  function,  r*,  is 
U-shaped.  The  value  of  .Ml  where  r*(..K)  is  minimal,  is  proposed  as  a 
quantitative  location  of  the  transition  from  slow  to  fast  growth, 
with  symbol  .\K 

i.e.,  for  . K - K , -Kj-  0.  (!>2) 


To  demonstrate  the  operation  of  the  previous  section's  analysis  an 
example  is  con-idered  using  T i -GAT - GV-ifbn  fatigue-crack  growth  data. 
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This  data  was  generated  and  reduced  at  University  of  llissouri- 
Columbia,  during  Winter  Semester  1975. 

5.3.2  numerical  Analysis  of  Crack  Propagation  Measurements 


5.3.2. 1 Purpose  of  Improved  Curve  Fitting 

Numerical  procedures  can  be  employed  for  predicting  the  number 
of  load  cycles  for  crack  growth  from  an  initial  crack  length  to  a 
final  length  corresponding  to  a design  life  or  an  inspection  interval, 
These  procedures  can  be  evaluated  using  compact  tension  (CT)  or  wedge 
open  loading  (WOL)  specimen  conpliance  calibration  and  crack  growth 
rate  test  results.  As  in  section  5.2.3  the  WOL  specimen  strpss  in- 
tensity factor  K can  be  expressed  by  rearranging  Lq.  31  as: 


y - Eli!,  f c (9.)  d ( CL  13 ) -il/2  _ i-P>  a-ip 
K " Ba  L V d(a?WTJ  ~ LBa  J 3 


(53) 


. - r c ,a\  d(CCB),l/2  . rM,  . -flB  . 

where  C^  r [.5  (^)  j-pp/jjy J and  Cl  B - y is 


the 


compliance  calibration  parameter.  Calculation  of  K requires  a dif- 
ferentiable fitting  relation  for  CCB  as  a function  of  a/W.  One 
approach  used  widely  in  the  past  by  fracture  mechanicians,  is  to 
apply  a compliance  calibration  expression  of  the  following  form: 


83 


6 i ^ 

CEB  = Exp[  l k.  (a/W)1’1]  (54) 

i = l 1 

usually  considered  valid  for  0.3  <_  a/W  <_  0.8.  The  polynomial  co- 
efficients kj  can  be  determined  wltn  the  aid  of  an  interactive  pro- 
gramming system  (PL/1)  curvilinear  regression  routine,  POLFIT. 

The  results  of  the  polynomial  regression  analysis  for  the 
calibration  data  can  be  used  in  conjunction  with  the  curve  fitting 
relation  to  express  the  polynomial  compliance  calibration.  This 
curve  segment  has  been  seen  to  fit  adequately  for  a range  in  a/W 
between  0.3  to  0.8.  Using  WOL  crack  propagation  results  as  a test 
case  life  or  inspection  interval  prediction  exercises  can  be  per- 
formed . 

It  Is  desirable  to  expand  the  boundary  conditions  so  that  the 
stress  Intensity  range  AK,  given  by 

AK  = (1  - load  ratio,  R)  C3,  (55) 

could  be  calculated  for  values  of  outside  the  range  0.3  to  0.8 
of  a valid  compliance  calibration  fitting  relation. 

Consider  a more  useful  fitting  relation  that  would  extend  the 
allowed  range  of  a/W  to  0 <_  a/W  <_  1.0.  To  explore  the  possibility, 
the  Wei  bull  distribution  was  determined  to  display  better  physical 
relation  to  the  required  boundary  conditions.  As  previously 
mentioned  an  interactive  (PL/1)  routine  called  WYBULL  was  developed 
(Bowie  et  al , 1974)  for  estimating  the  three  parameters  k,  e,  v of 
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the  Welbull  distribution  function: 

P(x)  = Exp  [-  (~~)k]  (56) 

Minor  modifications  of  the  routine  yielded  a new  program,  called 
WOLFIT  that  Is  used  to  obtain  the  three  parameters  k,  e,  and  v,  for 
the  new  WOL  specimen  compliance  calibration  fitting  relations  In 
section  5. 3. 1.2 

CEB  = H and  a/W  = 1 - F 

The  boundary  condition  of  a/W  = 1 

requires  1 - f * 0,  Is  satisfied  when  lim  (CEB)  = °»  . 

a/W  - 1 

The  boundary  condition  of  a/W  = 0 

requires  1 - F « 1,  Implies  logeloge(CEB  + Exp(l))  - e « 0 
or  CEB (a/'l  = 0)  = ExpExp(e)  - Exp(l). 

For  example,  If  Exp(e)  = 1,  that  Is  e = 0,  then  CEB  = 0 when  a/W  = 0. 
When  comparing  the  polynomial  form  of  C3  and  that  of  the  Welbull 
distribution  form,  it  was  seen  that  the  difference  of  Welbull  C3  was 
negligible  for  the  range  0.3  to  0.6  of  a/W,  but  was  20°  lower  at 
a/W  = 0.1  and  at  a/W  = 0.8  was  10%  lower  than  the  polynomial  C3 
(Figure  23). 

5. 3. 2. 2 Initial  Test  Data  Reduction 

The  compact  tension  specimen  test  results  studied  to  date  con- 
sist of  sets  of  crack  length  and  load  cycle  values  obtained  In  ex- 
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perlments  performed  at  a set  of  loads  with  sin'  ■ form,  where  load 
ratio,  R = 0.05.  Crack  length  measurements  were  made  on  both  sides  of 
the  specimen.  A typical  raw  data  set  studied  was  for  the  case  of 
specimen  thickness  B = 0.996  in.  (2.53  cm),  and  width  VJ  = 3.0  In. 

(7.62  cm).  Load  cycles  were  applied  to  a frequency  of  10  Hz,  and 
the  test  was  performed  at  room  temperature  in  laboratory  air  environ- 
ment. Note  that  Table  IV  contains  48  pairs  of  crack  length  measure- 
ments and  48  numbers  of  load  cycles  at  which  the  crack  length  mea- 
surements were  performed.  Experimental  values  of  crack  propagation 
rate  are  given  by 

* 10'6  |6? 

, ,0-6  ALLLlJJ.t  I.Zjl j - MM)  - ft«.2))/2 

i --  1,2,3  ....  (58) 

where  A(j,l),  A(j,2)  are  crack  lengths  measured  on  sides  1 and  2 of 
the  specimen  and  fl(j)  is  the  cumulative  number  of  load  cycles.  Ex- 
perimental values  of  stress  Intensity  range  are  given  by 

AK(1)  = (1-R)  ks1  *'1"  (MPa  i = l .2,3. . (59) 

Derivations  of  numerical  expressions  for  function  were  discussed 
in  section  5.2.3.  The  value  of  K>b  was  chosen  to  be  50  ksi  ,Tn  (54.9 
MPa  ,'fn). 

A summary  of  Initial  test  data  reduction  for  specimen  T-92B-1 
run  at  the  University  of  Missouri,  Winter  1975,  is  aiven  in  Table  V. 
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TABLE  IV 

FATIGUE-CRACK  GROWTH  RAW  DATA 


Crack  Length 

Side  1 

Side  2 

Max.  Load 

inch 

Inch 

Cycles,  N 

kips 

1 

1.0520 

1.0540 

0 

4.0 

2 

1.0581 

1.0589 

2000 

4.0 

3 

1.0622 

1.0666 

4000 

4.0 

4 

1.0694 

1.0720 

6000 

4.0 

5 

1.0745 

1.0815 

8000 

4.0 

6 

1.0830 

1.0870 

10000 

4.0 

7 

1.0886 

1.0906 

12000 

4.0 

8 

1.0905 

1.0925 

14000 

4.0 

9 

1.0980 

1.1004 

16000 

4.0 

10 

1.1025 

1.1055 

18000 

4.0 

11 

1.1153 

1.1193 

20000 

4.0 

12 

1.1250 

1.1350 

22000 

4.0 

13 

1.1395 

1 .1415 

24000 

4.0 

14 

1.1360 

1.1560 

26000 

4.0 

15 

1.1420 

1.1560 

28000 

4.0 

16 

1.1670 

1.1710 

30000 

4.0 

17 

1.1690 

1.1810 

32000 

4.0 

18 

1.1865 

1.1895 

34000 

4.0 

19 

1.1960 

1.1980 

36000 

4.0 

20 

1.2055 

1.2105 

38000 

4.0 

21 

1.2140 

1.2160 

40000 

4.0 

22 

1.2239 

1.2311 

42000 

4.0 

23 

1.2345 

1.2385 

44000 

4.0 

24 

1.2462 

1.2526 

46000 

4.0 

25 

1 . 2608 

1.2632 

48000 

4.0 

26 

1.2720 

1.2750 

50000 

4.0 

27 

1.2829 

1.2871 

52000 

4.0 

28 

1.2995 

1.2995 

54000 

4.0 

29 

1.3120 

1 . 3096 

56000 

4.0 

30 

1.3261 

1.3241 

58000 

4.0 

31 

1.3390 

1.3374 

60000 

4.0 

32 

1.3593 

1.3573 

62000 

4.0 

33 

1.3665 

1.3625 

64000 

4.0 

34 

1.3780 

1.3764 

66000 

4.0 

35 

1.3950 

1.3936 

68000 

4.0 

36 

1.4145 

1.4139 

70000 

4.0 
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TABLE  IV  (cont.) 


Crack  Length 

Side  1 

Side  2 

Inch 

Inch 

37 

1.4280 

1.4276 

38 

1.4645 

1.4635 

39 

1.4798 

1.4788 

40 

1.5831 

1.5823 

41 

1.6229 

1.6209 

42 

1.6598 

1.6588 

43 

1.6964 

1.6954 

44 

1.7323 

1.7303 

45 

1.7761 

1.7741 

46 

1.8295 

1 .8285 

47 

1.9011 

1.9011 

48 

2.0192 

2.0172 

Max.  Load 
Cycles,  N kips. 


72000 

4.0 

74000 

4.0 

76000 

4.0 

90000 

4.0 

92000 

4.0 

94000 

4.0 

96000 

4.0 

98000 

4.0 

100000 

4.0 

102000 

4.0 

1 04000 

4.0 

106000 

4.0 

88 


TABLE  V 

FATIGUE -CRACK  GROWTH  REDUCED  DATA 


4a 

AN 

Da/DN 

AK 

Max.  Load 

1 

0.0045 

2000 

2.25 

12.843 

4.0 

2 

0.0059 

2000 

2.95 

12.890 

4.0 

3 

0.0063 

2000 

3.15 

12.941 

4.00 

4 

0.0073 

2000 

3.65 

13.002 

4.00 

5 

0.0070 

2000 

3.50 

13.061 

4.00 

6 

0.0046 

2000 

2.30 

13.100 

4.00 

7 

0.0019 

2000 

0.95 

13.116 

4.00 

8 

0.0077 

2000 

3.85 

13.183 

4.00 

9 

0.0048 

2000 

2.40 

13.225 

4.00 

10 

0.0133 

2000 

6.65 

13.344 

4.00 

11 

0.0127 

2000 

6.35 

13.461 

4.00 

12 

0.0105 

2000 

5.25 

13.560 

4.00 

13 

0.0055 

2000 

2.75 

13.612 

4.00 

14 

0.0030 

2000 

1.50 

13.641 

4.00 

15 

0.0200 

2000 

10.00 

13.838 

4.00 

16 

0.0060 

2000 

3.00 

13.898 

4.00 

17 

0.0130 

2000 

6.50 

14.031 

4.00 

18 

0.0090 

2000 

4.50 

14.124 

4.00 

19 

0.0110 

2000 

5.50 

14.240 

4.00 

20 

0.0070 

2000 

3.50 

14.315 

4.00 

21 

0.0125 

2000 

6.25 

14.451 

4.00 

22 

0.0090 

2000 

4.50 

14.551 

4.00 

23 

0.0129 

2000 

6.45 

14.696 

4.00 

24 

0.0126 

2000 

6.30 

14.840 

4.00 

25 

0.0115 

2000 

5.75 

14.974 

4.00 

26 

0.0115 

2000 

5.75 

15.111 

4.00 

27 

0.0145 

2000 

7.25 

15.286 

4.00 

28 

0.0108 

2000 

5.40 

15.418 

4.00 

29 

0.0148 

2000 

7.40 

15.603 

4.00 

30 

0.0133 

2000 

6.65 

15.772 

4.00 

31 

0.0199 

2000 

9.95 

16.031 

4.00 

32 

0.0062 

2000 

3.10 

16.113 

4.00 

33 

0.0127 

2000 

6.35 

16.283 

4.00 

34 

0.0171 

2000 

8.55 

16.517 

4.00 

35 

0.0199 

2000 

9.95 

16.796 

4.00 

36 

0.0136 

2000 

6.80 

16.991 

4.00 

37 

0.0362 

2000 

18.10 

17.528 

4.00 

i 


4 

1 


» ■ 4 
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The  object  of  crack  propagation  rate  curve  fitting  Is  to  select  a 
relation  to  represent  the  experimental  results  in  the  Da/DN  and 
AK/K^  columns  of  Table  V. 

The  Type-A  exponential  crack  propagation  curve  fitting  rela- 
tions defined  in  section  5. 3. 1.2  were  applied  to  the  measurements 
of  Oa/ON  and  AK./K^  listed  in  Table  VI.  The  PL/I  routine  was  used  to 
estimate  the  Weibull  parameters  k,  e,  and  v.  The  regression  analysis 
results  are  summarized  in  Figure  24,  where  the  circled  points  have 
regression  variable  coordinates 

V • [-  loge(l  - AK/Kb)]1/k  (60) 

and 

X = b Da/ON  + a 

5. 3.2. 3 Type-A  Curve  Fitting  Relations  for  Ti-6Al-6V-2Sn,  Specimen 
T-92B-1,  Tested  in  45.9%  R.H.  Lab  Air  at  R = 0.05 

Application  of  Type-A  crack  propagation  curve  fitting  relatior 
was  made  for  the  case  of  a set  of  47  Da/DN  and  £K  measurements  for 
a Ti-6Al-6V-2Sn  WOL  specimen  tested  under  a 10  Hz  sine  wave,  with 
R = 0,05  in  ambient  laboratory  air  (R.H.  = 45r  , SO  F (26.7  C)).  The 
specimen  T-92B-1  dimensions  were  as  foV'ows  specimen  thickness, 

U = 0.996  in  (2.53  cm)  width,  K = 3.C  in  (7.62  cm)  and  initial  flaw 
(machined),  = 0.90  in  (2.29  cm). 


TABLE  VI 

MEASURED  AND  CALCULATED  AK 


THE  CRACK  GROWTH  PARAMETERS  K,  E,  V ARE 


0.55710334 


“3.97261906 


57.96824646 


MEASURED 

CALCULATED 

DADN 

DELTAK 

GELTAk 

2.25000 

12.84300 

12.13429 

2.95000 

12.89000 

12.77290 

3.15000 

12.94100 

12.94811 

3.65000 

13.00200 

13.37341 

3.50000 

13.06100 

13.24767 

2.30000 

13.10000 

12.18129 

0.95000 

13.11600 

10.82383 

3.85000 

13.18300 

13.53870 

2.40000 

13.22500 

12.27463 

6.65000 

13.34400 

15.61683 

6.35000 

13.46100 

15.41237 

5.25000 

13.56000 

14.62808 

2.75000 

13.61200 

12.59458 

1 . 50000 

13.64100 

11.40045 

10.00000 

13.83800 

17.67674 

3.00000 

13.89800 

12.81698 

6 . 50000 

14.03100 

15.51508 

4.50000 

14.12400 

14.05850 

5 . 50000 

14.24000 

14.81138 

3.50000 

14.31500 

13.24767 

6.25000 

14.45100 

15.34337 

4.50000 

14.55100 

14.05850 

6.45000 

14.69600 

15.48095 

6.30000 

14.84000 

15.37793 

5.75000 

14.97400 

14.99162 

5.75000 

15.11100 

14.99162 

7.25000 

15.28600 

16.01471 

5.40000 

15.41800 

14.73843 

7.40000 

15.60300 

16.11198 

6.65000 

15.77200 

15.61683 

9.95000 

16.03098 

17.64859 

3.10000 

16.11298 

12.90460 

6.35000 

16.28297 

15.41237 
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TABLE  VI  (cont.) 


8.55000 

16.51698 

16.83093 

9.95000 

16.79597 

17.64859 

6.80000 

16.99098 

15.71765 

18.09999 

17.52798 

21.51921 

7.65000 

17.76297 

16.27225 

7.39000 

19.50098 

16.10553 

19.59999 

20.23698 

22.11096 

18.70000 

20.98798 

21.75949 

18.29999 

21.77397 

21.59984 

17.70000 

22.58897 

21.35626 

21 .89999 

23.67998 

22.96436 

26.95000 

25.16798 

24.64584 

36.04999 

27.45598 

27.17189 

58.54999 

32.09697 

31 .70181 

II 


Figure  24  Regression  Analysis  for  Ti-6Al-6V-2Sn  Fatigue-Crack 

Growth  Data  for  a Mill  Annealed  High  Oxygen  Alloy  (Speci- 
men T-92B-1). 


i j 


■ i 

{ i 


i i 

1 

<2 

1 


1 


! 


* I i 
i i 


Li 


93 


94 


For  the  example,  the  value  of  Kb  was  estimated  to  be  50  ks 1 
An  (54.9  MPa  ym).  The  results  of  the  curve  fitting  exercise  are 
shown  In  Figure  25,  where  the  ordinate  is  Da/DN  to  a logarithmic 
scale  and  the  abscissa  Is  AK/Kfa  to  a linear  scale.  The  dash-dot 
curve  represents  the  median  fitting  relation 

Da/DN  = e + (v-e)(-loge(l  - AK/Kb))1/k  (61) 

where 

e = -3.9726 
v = 57.9682 
k = 0.5571 

Kb  = 50  ksi  An  (54.9  MPa  A\) 

The  threshold  stress  intensity,  Kth,  was  calculated  to  be  equal 
to  9.733  ksi  /Tn  (10.71  MPa  /m) . Measured  and  calculated  values  of 
AK  are  compered  In  Table  VI  at  the  47  values  of  measured  Da/DN.  The 
values  are  plotted  in  Figure  25  with  the  AK  being  normalized. 


5. 3. 2.4  Use  of  Risk  Analysis  to  Distingui:  w and  Fast  Crack  Growth 


For  the  Welbull  distribution  function,  used  in  fatigue  data 
analysis,  i.e. 

^ s>  Is 

(62) 


P(x)  = Exp  [-  • 


'v  - e 

It  Is  usual  to  define  a risk  function,  r*(x),  as 

_ k /X  - esk-1 
(v  - e)  - e; 


(63) 


CRACK  PROPAGATION  RATE,  oa/ON,  „ in/cycle 


lOOOl 


0.5571 

-3.9726 

57.9682 


10000 


/ -4  1000 


K,  ■ 52.5  ksl/fn 
xc  (57.6  MPa^S) 

K.  *>  (l-R)K,  - 50  ksl/Tn 

15  Ic  (54.9  MPayffi) 

Plotted  Sample  Size  N » 30 

Total  Sample  Correlation 
Coefficient,  r * 0.97242 


! Kth  ■ 9.733  ksl/Tn  (10.71  HPtJi) 


NORMALIZED  STRESS  INTENSITY  RANGE.  .WKb 


CRACK  PROPAGATION  RATE,  Ds/dn,  v cm/cyd 
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For  the  purposes  of  this  work,  a new  definition  of  risk  func- 
tion was  developed  which  holds  promise  of  being  useful  for  dis- 
tinguishing between  slow  and  fast  growth  segments  of  curve  fitting 
relations  A,  B and  C.  The  new  risk  function  Is  defined  as 


r* 


_] 6(Da/DN) 

Oa/DN  5(AK/Kb) 


For  example,  when  the  risk  function  and  slope 


(64) 


(65) 


are  known  at  some  AK/K^,  then 

WO"  ■ hr  <66> 

The  crack  propagation  curve  fitting  relation,  of  Type-A  was 
used  to  calculate  a numerical  risk  function,  summarized  in  Figure 
26.  Note  that  the  curve  is  U-shaped.  It  Is  proposed  that  the  value 
of  AK/Kb  at  the  minimal  values  of  r*  Is  a suitable  point  for  dis- 
tinguishing between  slow  and  fast  crack  growth  on  the  median  Oa/DN 
versus  AK/K^  curve. 


5.3.3  Numerical  Integration  of  Crack  Propagation  Relations 
5. 3. 3.1  The  Crack  Propagation  Integral 


ihe  main  reason  why  crack  propagation  test  results  are  ex- 
pressed l>y  relations  of  the  form 


Da/DN  f (AK) 


(67) 


Figure  26  Risk  Analysis  Plot  to  Determine  Slow  to  Fast  Transition 
in  Crack  Growth  Rate  Data. 
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Is  to  permit  prediction  of  the  number  of  cycles  for  growth  from  an 
initial  crack  length  a^  to  a final  length  af  consistent  with  a 
cyclic  load  history  that  Is  usually  different  than  the  load  his- 
tory for  which  experimental  crack  propagation  test  results  were 
obtained.  The  Increment  In  the  number  of  load  cycles  for  crack 
growth  from  length  a^  to  length  a^  Is  defined  by  the  Integral 

A 106 

MfNi*]  h7MrftKT  da>  cyc1eSt  (68) 

a1 

The  usefulness  of  a crack  propagation  rate  curve  fitting  relation 
can  be  judged  by  successive  Integration  using  the  sequences  of  com- 
pact tension  specimen  crack  lengths  as  limits  of  Integration,  while 
taking  due  account  of  the  test  load  history.  The  results  can  be 
plotted  In  the  form 

a = a(N) 

and  can  be  compared  with  experimental  pairs  (a,N). 

5. 3. 3. 2 Numerical  Integration  Routine 

The  numerical  form  of  the  crack  propagation  Integral  can  be 
calculated  with  the  aid  of  a PL/I  routine.  The  routine  employs  a 
variable  mesh  Simpson's  rule  Integration  scheme  (Bowie  et  al , 1974). 
The  argument  of  the  integral  is  first  computed  at  the  11  crack 
lengths, 
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a( j ) * a1  ♦ (j  - 1)  h.  In  (cm),  j = 1,2,. ...11,  (69) 

where 

h 8 («u  - a1 )/10, 

a^  = lower  limit  crack  length.  In  (cm) 
au  * upper  limit  crack  length.  In  (cm). 

At  each  crack  length  a(j),  the  normalized  stress  Intensity  range, 
AK/K^  Is  calculated  as  follows: 


The  values  of  Da/0N(j)  at  the  11  cracr.  lengths,  are  determined 
using  one  of  the  Type-A,  A,  o,'  C curve  fitting  relations  defined  in 
paragraph  5. 3. 1.2.  The  numerical  form  of  the  crack  propagation  In- 
tegral is  given  by 

iN  ' 10*tj  (BS7BHTTT  * fia/BHTu) } * J,2>18  T ffWjT 

* , r Da/BfTCJT3 * c),cles-  (71) 

The  mesh  size  Is  Increased  from  10  to  12  and  the  process  Is  re- 
peated at  13  crack  lengths 

a(j)  « a}  + h(j-l).  In  (cm),  j = 1,2,3,  ...  13  (72) 

where 


h = (ay  - a1 )/12 
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and  the  absolute  difference  between  AN  values  for  the  two  mesh 
sizes  Is  compared  with  a specified  error  Increment*  usually  taken 
to  be  1.0  cycle.  Iteration  of  the  numerical  procedure  Is  per- 
formed with  the  mesh  size  increased  by  two  until  the  specified 
precision  is  achieved. 

5. 3. 3. 3 Recovery  of  (a,N)  Data  by  Integration  of  Compact  Tension 
Specimen  Crack  Propagation  Curve  Fitting  Relation 

The  object  of  integrating  crack  propagation  relations  of  the 

form 

Da/DN  = f(AK) 

is  to  find  the  relation  between  crack  length,  a,  and  load  cycles, 

N,  for  a particular  load  history  to  predict  inspection  intervals 
and  estabMsh  operating  ranges.  The  condition  for  failure  is 
simply 

AK  = Kb  ksi  An  (MPa  At)  (73) 

For  the  exponential  crack  propagation  curve  fitting  relations, 
the  term 

[-  loge  (1  - AK/Kfc)]1/k  (74) 

is  singular  for  the  failure  condition.  In  practice,  when  applying 
the  numerical  integration  routine, the  failure,  or  unstable  crack 
growth,  condition  is  defined  as 


0.99  K.  < AK  < K.  ksi  /Tn  (MPa  r'm) 
D — D 
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It  would  seem,  in  theory,  possible  to  integrate  the  exponential 
Da/ON  relations  at  increasing  AK  beginning  at 

AK  * Kth  + c = Kb(l  - Exp  [-  (^-)]k)  + e,  ksl  /TF  (MPa  &)  (75) 

where  e is  arbitrarily  small.  When  this  approach  is  followed  at 
selected  e,  numerical  difficulties  are  often  encountered  when 
attempting  to  reproduce  measured  crack  length  and  load  cycle  data  sets. 
The  difficulties  relate  to  the  following  factors: 

1.  Resolution  of  crack  length  increments  at  AK 
values  near  Kth  Is  inferior  with  respect 

to  resolution  at  relatively  high  AK  values. 

2.  Crack  propagation  rate  experiments  at  low  AK 
are  relatively  time  consuming,  with  the  result 
that  median  crack  propagation  curve  fitting  re- 
lations are  cocrmonly  weighted  to  a minor  degree 
by  the  low  AK  measurements. 

3.  Measurements  at  low  AK  are  apt  to  reflect  local 
differences  In  test  specimen  microstructure,  with 
the  result  that  scatter  in  Da/DN  and  AK  measure- 
ments is  relatively  severe  and  Is  not  well  repre- 
sented by  median  curve  fitting  relations  determined 
by  analysis  of  the  entire  data  set  for  a specimen. 

At  present,  judgement  must  be  exercised  to  select  the  lowest  AK  at 
which  Integration  of  an  exponential  curve  fitting  relation  should 


be  attempted. 

The  results  of  applying  a PL/I  routine,  LIFE,  to  integration 
of  the  exponential  Oa/DN  relation  for  specimen  T-92B-1  with  two 
beginning  points  of  Integration  from  aQ  * 0.9  In  (2.29  cm)  and 
aQ  = 1.0  in  (2.54  cm)  are  plotted  In  Figure  27.  The  ordinate  re- 
presents load  cycles,  N,  to  a logarithmic  scale  and  the  abscissa 
indicates  the  crack  length,  a.  The  circled  points  are  experimental 
measured  values.  The  two  curves  are  numerical  Integration  results; 
Curve  1 assumes  a K associated  with  the  last  measureable  a^  = 2.0182 
in  (5.126  cm)  with  an  initial  crack  length  aQ  = 0.9  In  (2.29  cm)  and 
Curve  2 assumes  a equal  to  50  ksi  /TF  (54.9  MPa  v'm)  with  an  Initial 
crack  length  aQ  - 1.0  in  (2.54  cm)  and  a predicted  final  crack  length 
a^  = 2.1404  In  (5.437  cm).  The  choice  of  the  first  curve  and  close 
scrutiny  In  the  measurement  of  the  W0L  specimen  yield  the  most  accu- 
rate correlation.  The  measured  and  calculated  values  are  given  in 
Table  VII. 

5. 3. 3.4  Examples  of  Practical  Graph  Formats  for  Life  or  Crack 
Growth  Internal  Prediction 

For  all  of  the  example  curves,  the  initial  assumption  was  made 
that  load  cycle  counting  began  when  the  crack  length  was  equal  to 
1 .0  Inch  (2.54  cm). 

Constant  Load  a-N-P  Curves  - Figure  28  Includes  8 curves,  for 
constant  load  amplitudes  from  P = 3 kips  (13.344  nt)  to  P = 10  kips 


Figure  27  Curve  Fit  for  Ti-6Al-6V-2Sn  Fatigue-Crack  Length  Versus 
Load  Cycles  Data  (Specimen  T-92B-1).  Curve  1 
illustrates  comparison  of  an  integration  scheme  using 
Simpson's  mesh  with  Wei  bull  parameters  and  measured  data 
for  aQ  = 0.9  inch.  Curve  2 plots  the  measured  data 

against  a predictive  technique  using  the  Weibull  para- 
meters, k,  e,  v,  and  Kb  = 50  with  an  initial  crack 

length  of  aQ  = 1.0  inch. 
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CRACK  LENGTH,  a,  cm 


0.9  1.1  1.3  1.5  1.7  1.9  2.1 

CRACK  LENGTH,  a,  inch 
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TABLE  VII 

MEASURED  AND  CALCULATED  FATIGUE  CYCLES 


Measured  Quantities 

Calculated 

Crack  Length,  In 

Cycles 

Cycles 

1 

.900 

' 0 

0 

2 

.904 

2000 

1495 

3 

.910 

4000 

3426 

4 

.917 

6000 

5450 

5 

.924 

8000 

7748 

6 

.931 

10000 

9904 

7 

.936 

12000 

11297 

8 

.937 

14000 

11867 

9 

.945 

16000 

14142 

10 

.950 

18000 

15534 

11 

.963 

20000 

19288 

12 

.976 

22000 

22735 

13 

.986 

24000 

25487 

14 

.992 

26000 

26895 

15 

.995 

28000 

27652 

16 

1.015 

30000 

32536 

17 

1.021 

32000 

33946 

18 

1.034 

36000 

36914 

19 

1.043 

36000 

38904 

20 

1 .054 

38000 

41265 

21 

1.061 

40000 

42728 

22 

1.073 

42000 

45268 

23 

1.082 

44000 

47040 

24 

1.095 

46000 

49500 

25 

1.108 

48000 

51817 

26 

1.119 

50000 

53859 

27 

1.131 

52000 

55834 

28 

1.145 

54000 

58235 

29 

1.156 

56000 

59960 

30 

1.171 

58000 

62239 

31 

1.184 

60000 

54208 

32 

1.204 

62000 

67021 

33 

1.210 

64000 

67866 

34 

1.233 

66000 

69552 

35 

1.240 

68000 

71730 

36 

1.260 

70000 

74139 

37 

1.274 

72000 

75711 

38 

1.310 

7400C 

79619 
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TABLE  VII  (cont.) 


39 

1.325 

76000 

40 

1.429 

90000 

41 

1.468 

92000 

42 

1 .505 

94000 

43 

1.542 

96000 

44 

1.577 

98000 

45 

1.621 

100000 

46 

1.675 

102000 

47 

1.747 

104000 

48 

1.864 

106000 

Measured  Quantities 
Crack  Length,  In  Cycles 


1 

1.0 

0 

2 

1.015 

4884 

3 

1.021 

6294 

4 

1.043 

11252 

5 

1.050 

6 

1.061 

15076 

7 

1.082 

19388 

8 

1.100 

9 

1.108 

24164 

10 

1.119 

26204 

11 

1.145 

30583 

12 

1.156 

32308 

13 

1.184 

36556 

14 

1.200 

15 

1.204 

39369 

16 

1.210 

40214 

17 

1.240 

44078 

18 

1.260 

46487 

19 

1.300 

20 

1.310 

51967 

21 

1.400 

22 

1.429 

62367 

23 

1.468 

65113 

24 

1.500 

25 

1.505 

67462 

26 

1.542 

69525 

81158 

90018 

92765 

95114 

97176 

98968 

100929 

102990 

105209 

107714 


Calculated 

Cycles 

0 


11514 


21434 


37525 


.49675 

58904 


65892 


16 


77 

38 

57 

61 

79 


74923 

77590 

79223 

79952 

80713 

80909 

80956 

80981 

80982 


CRACK  LENGTH,  a,  cm 
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Initial  Crack  length  a.  » 1.0  In 
1 (2.54  cm) 

Stress  Ratio  R ■ 0.05 
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(44,480  nt)  at  1.0  kip  (4,448  nt)  Interve'*  The  locus  labelled 

> 

Da/DN  = 69.8  u In/cycle  (1772.9  u mm/cycle)  separates  slow  and  fast 
growth  regimes,  determined  by  the  risk  function  analysis  method  of 
section  5.3.1 .5. 

Constant  Crack  Length  N-P-a  Curves  - Figure  29  Includes  4 
curves,  for  constant  crack  length  from  a * 1.1  In  (2.79  cm)  to 
a = 1.4  In  (3.56  cm)  at  0.1  in  (0.254  cm)  Intervals.  The  unstable 
crack  growth  locus  Is  reminiscent  of  a conventional  low  cycle  fati- 
gue S-N  curve.  The  curve  separating  slow  and  fast  growth  regimes 
was  determined  by  the  risk  analysis  method. 

Load  Shedding  Curves  - The  staircase  function  plotted  In  Figure 
29  represents  a load  shedding  example  In  which  load  cycle  counting 
began  at  crack  length  a = 1.0  Inch  (2.54  cm)  with  load  P = 10  kips 
(44,480  nt).  The  load  was  reduced  by  1.0  kip  (4,448  nt)  Increment 
for  each  0.1  Inch  (0.254  cm)  Increment  in  crack  growth  until  a repre- 
sentative lower  load  limit  P = 3 kips  (13,344  nt)  was  attained.  This 
lower  load  limit  was  maintained  in  the  numerical  analysis  until  pre- 
dicted unstable  crack  growth  was  attained.  The  same  example  Is  re- 
presented In  Figure  28  by  the  curve  labelled  "A". 

The  curve  segment  labelled  "Aa"  In  Figure  28  was  generated  to 
Illustrate  an  Important  principle  for  In-service  structural  component 
life  prediction  when  load  shedding  Is  practiced.  The  numerical  analy- 
sis for  the  above  load  shedding  example  was  repeated  until  crack 
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nt).  The  load  was  then  Increased  to  P * 7 kips  (31,136  nt)  for  growth 
from  a 8 1.5  Inch  (3.81  cm)  to  a 8 1.55  Inch  (3.94  cm).  The  load 
was  then  further  Increased  to  P s 8 kips  (35,584  nt)  and  was  main- 
tained until  predicted  unstable  growth  was  attained.  This  exercise 
Illustrates  the  principle  that  when  load  shedding  Is  practiced  for 
an  In-service  structural  component,  a decision  to  Interrupt  a load 
shedding  schedule  In  favor  of  brief  application  of  higher  than 
scheduled  loads  can  result  in  drastic  reduction  In  component  life. 

Two  additional  load  shedding  examples  are  sunmarlzed  In  Figure 
30.  Comparison  of  curves  B and  C reveals  the  effects  of  delaying 
load  shedding  on  crack  growth  and  life  prediction.  Let  curve  B re- 
present an  in-service  structural  component  nominal  load  shedding  ex- 
ample, and  curve  C a possible  deviation.  If  the  two  components  were  \ 

4 

1 

to  be  Inspected  after  7000  load  cycles  had  been  applied,  predicted  \ 

crack  growth  for  curve  B would  be  from  0.9  Inch  (2.29  cm)  to  1.252 

Inch  (3.18  cm)  and  from  0.9  Inch  (2.29  cm)  to  1.500  Inch  (3.81  cm)  j 

3 

1 

for  curve  C.  If  the  next  Inspection  were  to  be  delayed  until  10,000  \ 

load  cycles  had  been  applied,  the  crack  length  for  curve  B would  be 

1.404  Inch  (3.57  cm),  and  1.72  Inch  (4.37  cm)  for  curve  C.  At  this  < 

stage,  component  B would  have  a remaining  life  of  approximately 

15,000  load  cycles,  and  component  C approximately  6,000  cycles. 
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Figure  30  Loading  Example  to  Demonstrate  Component  Shedding  Regimes 
as  Applied  to  Specimen  T-92B-1.  Curve  "B"  illustrates 
a load  P = 10  kips  for  973.7  cycles,  with  slow  growth 
from  10  kips  to  3 kips  determined  by  a scheme  that  pro- 
hibits fast  growth  rates.  From  N = 973  to  23821  cycles 
the  risk  function,  r*,  has  a minimum  value,  r*  = 

5.1823.  For  23821 <N<25375  cycles,  fast  growth  occurs 
at  load  P = 3 kips.  For  curve  "C",  the  initial  load 
P - 10  kips  is  maintained  until  a = 1.163  inches, 
at  N = 3205  cycles.  Load  shedding  then  takes  place, 
from  P * 10  kips  to  P = 3 kips  with  an  average 
minimum  risk  value,  r*  = 5.552. 
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CHAPTER  6 


TECHNICAL  DISCUSSION  OH  STUDY  RESULTS 

A graphical  sequence  of  analysis  in  fracture  and  fatigue  life 
mechanics  was  explored  with  respect  to  a single  alloy  test  in 
chapter  5.0.  Recognizing  the  discrimination  potential,  this 
study's  major  thrust  was  to  evaluate  differences  in  materials' 
microstructure  through  their  response  to  various  combinations  of 
parameters  as  applied  through  fatigue-crack  propagation.  The  ap- 
proach was  to  perform  fatigue-crack  propagation  testing  on  standard 
WOL  specimens,  at  three  different  laboratories,  and  use  a "best- 
fit"  approach  to  data  as  determined  by  the  definition  of  mathema- 
tical boundary  conditions  to  the  physics  of  fatigue-crack  propaga- 
tion (section  5. 3. 1.2).  It  appears  that  this  idea  coupled  with  sub- 
stantial and  precise  data  collection  for  quantizing  factors  such  as 
compliance,  crack  length  and  K,  matched  with  this  data  reduction 
technique,  resulted  in  distinguishable  data  trends  for  heat  treat- 
ments, stress  ratios,  wave  forms  and  environments  applied  to  the 
test  materials.  The  results  are  represented  by  two  fatigue-crack 
growth  rate  (DaDH)  versus  stress  intensity  (AK)  plots,  one  log- 
linear,  the  other  log-log.  Both  plots  reflect  the  Weibull  charac- 
teristic value  v,  intersecting  the  log  ordinate  axis  (DaDN)  and  a 
value  equivalent  to  0.63212  times  the  final  «b,  intersecting  the 
linear  abscissa.  All  comparative  analyses  were  performed  with  log- 
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linear  representations.  The  latter  curves  (Appendix  1)  were  In- 
cluded for  traditional  benefit  and  to  provide  a bridge  from 
straight  line  expressions  to  a suitable  curve  fitting  formula- 
tions that  attempt  to  establish  boundary  conditions  In  accordance 
with  the  mathematical  physics  of  the  event.  It  can  be  noted  that 
all  graphs  have  curves  which  asymptotically  approach  a final 
value  Kfa.  Conversely,  a threshold  can  be  considered  as  a 
point  where  the  fatigue-crack  growth  rate  DaDN  is  low  enough  to  be 
considered  virtually  unidentifiable  at  the  present.  Equation  41 
should  not  be  misconstrued  as  DaDN  = 0 even  though  It  uses  this 
method  to  arrive  at  a limiting  Kth>  Alternatively,  it  is  best  to 
be  considered  as  a value  for  DaDN  approaching  zero. 

Another  graphical  format  similar  to  Figure  25  plots  the  Welbull 
analysis  reduced  to  one  shaping  parameter,  k.  Here  the  log  ordinate 
Is  the  Z variate,  (X-e)/(v-e)  and  the  linear  abscissa  Is  the  norma- 
lized stress  Intensity  range,  AK/Kb> 

6.1  COMPLIANCE  EVALUATION 

The  Wilson  polynomial  (section  5.2.3)  was  used  in  the  calibra- 
tion of  K permitting  comparison  with  future  data.  Albeit  this  ex- 
pression does  not  attempt  to  approximate  an  ever  changing  Modulus  of 
Elasticity,  nor  plastic  zone  size,  but  at  the  present  It  does  en- 
deavor to  provide  continuity  to  a specimen's  stress  intensity  factor, 
K.  All  base-line  specimens  were  subject  to  compliance  testing.  Cal- 
culations using  an  optically  judged  crack  length  (section  4.4)  (4.e. 
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natural  fatigue  crack  versus  a machined  notch)  were  made  and  plotted 
In  Figure  31.  The  data  aligned  between  test  specimens  and  was  fur- 
ther conpared  with  data  available  from  Sullivan  (1975).  This  latter 
comparison  presented  in  Table  VIII  of  Figure  31  uses  machined  cut 
notches.  Agreement  is  seen  at  the  lower  a/W  ratio  and  deviated 
slightly  at  higher  a/W  ratios.  This  may  be  due  to  subtle  crack 
front  eccentricity  or  the  use  of  natural  fatigue  crack. 

6.2  WE I BULL  PARAMETER 

The  logic  for  estimating  the  three  parameters  (k,e,v)  is  with 
an  interactive  program  (section  5. 3. 1.2).  The  parameters  are  derived 
from  the  linearization  of  the  Weibull  equation  (equation  38'1  de- 
ducing the  appropriate  values  of  W,  a,  and  b.  Optimization 
data  fit  by  linear  regression  analysis  results  by  a numeral  anal- 
ysis "best-fit"  when  a maximum  sample  correlation  coefficient  is 
reached.  The  rationale  for  implementing  a Weibull  type  fitting  tech- 
nique lies  in  the  Weibull 's  power  to  segregate  subtle  changes  in 
data  sets.  This  capacity  to  describe  behavior  by  only  three  and 
possibly  a single  parameter  makes  it  attractive.  Depending  on  an 
alloy's  microstructure,  what  is  considered  a discreet  selection  or 
change  in  operational  variables,  based  on  existing  data  basis  can 
induce  substantial  unsuspected  limitations  or  lead  to  total  dete- 
rioration of  the  structure.  Figures  29  and  30  are  indicative  of  a 
typical  scenario  that  would  ignore  a load  shedding  regime.  Identi- 
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Figure  31  Compliance  Evaluation  Curve.  The  compliance  response  used 
for  this  study  Is  plotted  with  the  data.  Comparison  Is 
shown  In  the  Incorporated  Table  VIII  that  presents  CEB 
data  for  the  University  of  Missouri  and  Naval  Research 
Laboratory  at  the  Indicated  a/W  levels. 
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fication  of  the  responses  to  adjustments  in  operation  variables  may 
be  recognized,  but  the  differences  or  deviations  may  appear  muted 
in  rapid  or  limited  laboratory  tests  and  totally  non-existent  when 
results  are  plotted  in  a log-log  regime.  Thus,  a needed  capability 
has  been  recognized  not  only  in  an  academic  sense  but  has  gained 
visibility  by  the  continuing  limitations  placed  on  field  operation 
of  high  performance  structures. 

6.2.1  Fitting  Data  for  Test  Matrix 

A fractional  factorial  test  matrix  was  pursued  for  both  ma- 
terials per  Table  IX.  The  resulting  fitting  data  associated  with 
these  data  sets  are  shown  in  the  bottom  portion  of  Table  IX.  As 
will  be  seen  the  fatigue-crack  growth  rate  curves  are  plotted  with 
a characteristic  value,  v,  against  a normalized  abscissa  AK/K^. 

This  can  be  interpreted  as  the  transition  point  betv/een  slow  and 
fast  crack  growth.  Further,  the  characteristic  can  be  used  as  a 
standardized  point  or  fulcrum  in  the  following  plots  which  look  at 
the  shaping  factor. 

By  defining  Z = (X-e)/(v-e)  as  a reduced  variate,  where  X = 

1/ 

DaDtl,  the  survivorship  function,  P(Z)  - Exp[-Z  ],  can  be  related 
to  the  probability  of  failure  by  Equation  35.  Thus,  l-P(Z)  Is 
identified  as  1-F  or  1-AK/Kb;  hence,  P(Z)  = 1-Exp[-Zk].  It  can  be 
noted  if  Z = 1,  AX/K^  = l-Exp(-l)  = 0.63212  and  when  1 = 0,  then 
AK/Kb  = 0 for  all  k.  Using  7 and  AK/Xb  as  axes  the  "Spider"  plots 
in  Figures  32  and  33  were  produced  to  illustrate  the  generic  be- 
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TABLE  IX  WEIBULL  TYPE-A  FITTING  PARAMETERS 


Sine  Wave,  10  Hz 
Ambient  Air 


0.05 


1.  T* 
8,  H* 


1 min  Hold-Time 
Ambient  Air 


0.05 


2 

12,  N* 


0.67 


6,  E* 
K* 


1 min  Hold-Time 
Sea  Water 


10,  B* 


7-4  PH 


1 

1.8518 

-23.3152 

51.0053 

0.0375 

0.9780 

16.5  150 

2 

1.4285 

-194.8995 

218.3016 

0.1873 

0.7710 

49.5  170 

5 

0.5000 

-1.4917 

16.4075 

0.1310 

0.9118 

12.5  50 

6 

0.5310 

-2.8944 

15.7506 

0.9292 

0.8004 

18.9  61 

8 

1.2853 

-17.2463 

57.8724 

0.0195 

0.9948 

21.0  150 

10 

0.4273 

-0.1453 

177.0029 

0.1440 

0.9540 

7.5  160 

11 

0.5000 

-0.8699 

15.8099 

0.1108 

0.9711 

12.3  60 

12 

15-5  PH* 

0.9352 

-21.9354 

100.1449 

0.0736 

0.9865 

30.9  170 

.8368 

.7813 

.6083 


0.6968 


i.534 

.923 


24.4371 

16.2135 

-6.5903 

65.4556 


0.2630 


.3461 

.2386 


.9821 
,7828  27.9  63 
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Figure  32  A "Spider"  Diagram  of  17-4  PH  for  all  Test  Conditions  In 
this  Study  Is  Graphed.  The  reduced  variate,  Z * (DaDN-e)/ 
(v-e)  Is  plotted  against  the  normalized  stress  Intensity 
AK/Kb  to  observe  the  shaping  factor  (K)  effect. 
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Figure  33  A "Spider"  Diagram  of  15-5  PH  for  all  Test  Conditions 
In  this  Study  Is  Graphed.  The  reduced  variate, 

Z * (DaDN-e)/(v-e)  Is  plotted  against  the  normalized 
stress  Intensity  AK/K.  to  observe  the  shaping  factor 
(K)  effect.  ° 
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havior  of  k. 

6.2.2  Behavior  Evaluation  by  Shaping  Factor 

In  considering  the  Influence  of  curve  shaping,  the  differences 
between  test  data  sets  must  be  Identified  and  partitioned  with  re- 
spect to  the  conditions  of  interest.  In  deriving  a rational  expla- 
nation for  overall  data  response  as  associated  with  a fitting  para- 
meter k,  and  in  the  attempt  to  align  a hypothesis  for  optimized 
fatigue-crack  growth  behavior,  five  distinct  sets  were  reviewed. 

These  data  :ets  that  are  Illustrated  and  compared  on  Figures  32 
and  33  are  1 and  8,  2 and  12,  5 and  11,  E and  K,  and  T and  H.  By 
observation  it  was  determined  that  to  optimize  (l.e.  achieve  maxi- 
mum life)  an  alloy's  fatigue  crack  behavior  two  items  are  required. 
First,  k must  tend  to  increase  positively,  secondly,  e must  tend  to 
increase  negatively. 

The  assumptions  considered  are:  that  initial  manufacturing  or 

service  flaws  can  be  Identified  for  the  materials'  forms  of  interest; 
and  the  test  conditions  reflect  the  actual  design's  operational  per- 
formance. With  these  boundary  conditions  it  becomes  essential  In  a 
llflng  analysis  to  Insure  that  entry  into  the  crack  growth  rate 
curves  does  not  occur  by  selective  internal  microstructural  bursting. 
This  last  factor  has  the  ability  to  erradicate  all  safety  margins 
along  with  some  equipment  and  structures. 

As  k increases  it  appears  to  indicate  a higher  resistance  to 
fatigue  damage  by  extending  the  curve  to  a reduced  crack  growth  rate 


over  a larger  stress  intensity  range.  This  can  be  seen  in  Figures 
32  and  33  for  Specimens  1 and  T,  respectively.  However,  e,  the 
threshold  value,  must  tend  to  be  more  negative  in  order  to  increase 
^th*  combination  is  a first  step  in  material  optimization  by 

data  set  scanning.  The  obvious  limitation  is  in  the  ability  to  define 
all  conditions  of  interest  and  compare  these  conditions  to  the  avail- 
able collected  data  set  of  interest.  This  study  indicates  the  sub- 
stantial impact  a set  of  applied  variables  have  on  alloys'  perfor- 
mance. 

Lastly,  the  characteristic,  v,  is  normalized  in  the  "Spider" 
diagrams  (Figures  32  and  33)  and  tends  to  act  as  the  pivot  point 
discriminating  between  slow  and  fast  crack  growth  rates.  Without 
definition  of  k,  e and  the  final  Kb,  the  interpretation  of  v is 
nebulus,  at  best.  However,  if  the  set  was  defined,  v has  the  poten- 
tial to  link  directly  with  design  methodology  safety  margins  for 
finite  life  designs. 

6.2.3  Collected  Uata  with  Curve-Fit 

It  is  believed  worthy  to  illustrate  in  log-linear  format  all 
the  fatigue-crack  growth  curves  of  the  data  set  with  their  data. 

The  Intercept  of  the  curve  indicates  a characteristic  value,  v, 
selected  to  be  the  slow  to  fast  transition  point  partitioning  the 
curve  into  two  distinct  parts.  The  curves  are  presented  in  numeric 
order  for  17-4  PH  and  alphabetic  order  for  15-5  PH  in  the  same  man- 
ner as  Table  IX  (Figures  34  to  48).  Identical  data  plots  made 


Figure  34  Log-Linear  Graph  of  Specimen  1 Fatigue-Crack  Growth  Rate 
Test  Fitted  with  a Type-A  Welbull  Curve.  The  slow  to  fast 
transition  Intercept  on  the  curve  indicates  a DaDN 
characteristic  value  (51.0  In/cycle)  and  63.21%  of  the 
final  stress  Intensity,  (150  ksi  /fn). 
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Figure  35  Log-Linear  Graph  of  Specimen  2 Fatigue-Crack  Growth  Rate 
Test  Fitted  with  a Type-A  Welbull  Curve.  The  slow  to 
fast  transition  Intercept  on  the  curve  Indicates  a DaDN 
characteristic  value  {218.3  In/cycle)  and  63.21%  of  the 
final  stress  intensity,  (170  ksl  /TrT). 


Figure  36  Log-Linear  Graph  of  Specimen  5 Fatigue-Crack  Growth  Rate 
Test  Fitted  with  a Type-A  Welbull  Curve.  The  slow  to 
fast  transition  Intercept  on  the  curve  Indicates  a DaDN 
characteristic  value  (16.4  In/cycle)  and  63.2)%  of  the 
final  stress  Intensity > (50  ksl  /fn). 
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Figure  37  Log-Linear  Graph  of  Specimen  6 Fatigue-Crack  Growth  Rate 
Test  Fitted  with  a Type-A  Wei  bull  Curve.  The  slow  to 
fast  transition  Inter rept  on  the  curve  indicates  a DaDN 
characteristic  val  ie  (15.7  In/cycle)  and  63. 2H  of  the 
final  stress  Intensity,  (61  ksl  /Tn). 
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Figure  38  Log-Linear  Graph  of  Specimen  8 Fatigue-Crack  Growth  Rate 
Test  Fitted  with  a Type-A  Melbull  Curve.  The  slow  to 
fast  transition  Intercept  on  the  curve  Indicates  a DaDN 
characteristic  value  (F7.8  In/cycle)  and  63.21%  of  the 
final  stress  Intensity,  (150  ksl  /Tn). 
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Figure  39  Log-Linear  Graph  of  Specimen  10  Fatigue-Crack  Growth  Rate 
Test  Fitted  with  a Type-A  Wei bull  Curve.  The  slow  to 
fast  transition  Intercept  on  the  curve  Indicates  a DaDN 
characteristic  value  (177.0  in/cycle)  and  63. ZH  of  the 
final  stress  Intensity,  (160  ksl  /Tn). 
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Figure  40  Log-Linear  Graph  of  Specimen  11  Fatigue-Crack  Growth  Rate 
Test  Fitted  with  a Type-A  Welbull  Curve.  The  slow  to 
fast  transition  Intercept  on  the  curve  Indicates  a DaDN 
characteristic  value  {15.8  In/cycle)  and  63.21%  of  the 
final  stress  Intensity,  {60  ksl  /Tn). 
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Figure  41  Log-Linear  Graph  of  Specimen  12  Fatigue-Crack  Growth  Rate 
Test  Fitted  with  a Type-A  Welbull  Curve.  The  slow  to 
fast  transition  Intercept  on  the  curve  Indicates  a DaDN 
characteristic  value  (100.1  In/cycle)  and  63.21%  of  the 
final  stress  Intensity,  1^  (170  ksl  /Tn). 
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Figure  42  Log-Linear  Graph  of  Specimen  B Fatigue-Crack  Growth  Rate 
Test  Fitted  with  a Type-A  Welbull  Curve.  The  slow  to 
fast  transition  Intercept  on  the  curve  Indicates  a DaDN 
characteristic  value  (221.8  In/cycle)  and  63.21%  of  the 
final  stress  in‘"is1ty,  Kb  (186  ksl  /Fn). 
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£ Figure  43  Log-Linear  Graph  of  Specimen  E Fatigue-Crack  Growth  Rate 

£ Test  Ftted  with  a Type-A  Welbull  Curve.  The  slow  to 

f fast  transition  Intercept  on  the  curve  indicates  a DaDN 

X characteristic  value  (15.9  In/cycle)  and  63.21%  of  the 

£ final  stress  intensity,  (63  ksi  /fn). 


Figure  44  Log-Linear  Graph  of  Specimen  H Fatigue-Crack  Growth  Rate 
Test  Fitted  with  a Type-A  Welbull  Curve.  The  slow  to 
fast  transition  Intercept  on  the  curve  indicates  a DaON 
characteristic  value  (151.6  In/cycle)  and  63.21%  of  the 
final  stress  Intensity,  (179  ksl  /Tn). 
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Figure  45  Log-Linear  Graph  of  Specimen  K Fatigue-Crack  Growth  Rate 
Test  Fitted  with  a Type-A  Welbull  Curve.  The  slow  to 
fast  transition  Intercept  on  the  curve  Indicates  a DaDN 
characteristic  value  (18.5  In/cycle)  and  63.21%  of  the 
final  stress  Intensity,  (60  ksl  */irT) . 
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Log-Linear  Graph  of  Specimen  L Fatigue-Crack  Growth  Rate 
Test  Fitted  with  a Type-A  Welbull  Curve.  The  slow  to 
fast  transition  Intercept  on  the  curve  Indicates  a DaDN 
characteristic  Vulue  (43.2  In/cycle)  and  63.21%  of  the 
final  stress  Intensity,  (58  ksl  /TrT). 
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Figure  47  Log-Linear  Graph  of  Specimen  N Fatigue-Crack  Growth  Rate 
Test  Fitted  with  a Type-A  Welbull  Curve.  The  slow  to 
fast  transition  Intercept  on  the  curve  Indicates  a DaDN 
characteristic  value  (163.9  In/cycle)  and  63.21%  of  the 
final  stress  Intensity,  (165  ksl  /Tn). 


Figure  48  Log-Linear  Graph  of  Specimen  T Fatigue-Crack  Growth  Rate 
Test  Fitted  with  a Type-A  Wei  bull  Curve.  The  slow  to 
fast  transition  Intercept  on  the  curve  indicates  a DaDN 
characteristic  value  (46.2  in/cycle)  and  63. 2U  of  the 
final  stress  Intensity,  Kfa  (150  ksi  /in). 
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straight  are  found  in  log-log  form.  Appendix  1. 

6.3  COMPARATIVE  ANALYSIS  OF  FATIGUE-CRACK  GROWTH  RATES 

The  major  point  to  be  derived  in  this  study  is  to  show  that 
micrcstructure  does  effect  fatigue-crack  growth  and  that  the  effect 
can  be  ever  entire  rate  and  stress  intensity  ranges.  The  primary 
thrust  is  to  examine  a result  that  may  have  its  explanation  and 
final  control  in  understanding  a combination  of  metal  physics, 
micromechanics,  continuum  and  fracture  mechanics.  Microstructural 
morphology,  tessellated  residual  strains,  secondary  particle  com- 
patibility and  activation  mechanisms  which  are  slightly  understood, 
do  direct  the  rate  at  which  materials  shall  be  reclaimed  to  a 
lower  energy  state  in  nature. 

To  appreciate  the  differences  between  the  following  test  set 
parameters  and  facilitate  viewing,  composite  graphs  were  con- 
structed. 

6.3.1  iieat  Treatment  Effect 

The  alloys'  microstructures  were  alterred  by  the  final  temper- 
ing or  aging  portion  of  the  heat  treat  process.  Both  alloys  had 
their  microstructures  put  in  the  H 1060  and  .1  1100  heat  treated  con- 
dition. The  basic  chemical  differences  (Table  III)  left  17-4  PH 
with  some  delta  ferrite  that  formed  into  extended  pools  and  stringers, 
whereas  16-5  PH  is  virtually  delta  ferrite  free.  The  retained 
austenite  was  too  fine  for  identification  in  either  alloy  or  heat 


treated  condition. 
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Initially  for  a 17-4  PH  base-line,  both  heat  treats  were  tested 
under  a 10  Hz  sine  wave  for  two  R-ratios,  R = 0.05  and  R = 0.67. 

Table  IX  will  identify  the  location  of  all  test  specimens  with  re- 
spect to  the  test  conditions  and  curves.  In  Figure  49,  the  test  ma- 
terial with  R = 0.05  and  H 1050  (curve  1)  has  a lower  threshold  than 
curve  8 at  H 1100.  But  as  curve  1 moves  toward  its  characteristic,  v 
its  crack  growth  rate  decreases.  For  R = 0.67  (curves  5 and  11), 
the  differences  were  negligible  for  threshold,  but  the  H 1100 
(curve  11)  tended  to  require  more  stress  intensity  for  fatigue 
cracking.  This  probably  is  not  surprising  in  that  the  H 1100  is  more 
tempered  with  a lower  yield  strength  and  a higher  resistance  to 
fracture.  However,  the  opposite  response  occurs  for  curves  1 and 
8.  In  Figure  50,  the  15-5  H 1050  and  H 1100  curves  T and  H align 
similarly  to  the  corresponding  17-4  curves  1 and  8 whereas  for 
R = 0.67,  the  15-5  H 1050  and  H 1100  curves  E and  K are  opposite 
to  the  17-4  curves  5 and  11,  and  the  15-5  curves  T and  H.  Further 
testing  at  R = 0.05  under  a one  minute  hold-time  trapazoidal  wave 
showed  distinct  mlcrostructural  response  differences  from  the  above 
mentioned  base-line.  At  H 1050  (curve  2),  the  threshold  In  conjunc- 
tion with  overall  fatigue-crack  growth  rates  increased  substantially, 
for  the  H 1100  (c^rve  12)  material  under  the  same  regime  neither 
threshold  or  rate  increased  as  much  as  curve  2.  However,  curve  12, 
still  exceeded  the  equivalent  10  Hz  base-line  of  curve  8.  It  can  be 
concluded  that  the  effect  of  heat  treatment  shown  by  the  microstruc- 
ture  can  vary  with  mean  stress.  Note  that  on  comparing  only  those 
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Figure  50  Composite  Log-Linear  Plot  of  Fatigue-Crack  Growth  Rates  of 
15-5  PH  Used  to  Relate  Heat  Treatment  (H  vs.  T,  K vs.  E). 
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test  variables  under  10  Hz  sine  wave,  the  15-5  PH  alloy  does  indi- 
cate slower  fatigue-crack  growth  rates  than  the  17-4  PH  alloy. 

6.3.2  R-Ratlo  Effect 

As  the  mean  stress  increases  the  fatigue-crack  growth  rates 
increase.  The  extent  of  the  increase,  is  largely  determined  on  the 
applied  stress  the  microstructure  experiences  (i.e.  including  its 
internal  residual  stress).  A possible  explanation  for  17-4  PH  H 1050 
(curve  5)  requiring  less  stress  intensity  to  fracture  may  be  found 
in  the  metal's  microscopy  and  fractography.  A more  precipitated 
matrix  with  larger  delta  ferrite  stringers  are  indicated  in  H 1100 
versus  H 1050.  This  combination  appears  to  be  an  active  retardant 
to  fatigue-crack  growth  over  a limited  AK  range  in  17-4  PH.  This  is 
seen  in  the  upper  portion  of  (curve  11),  and  in  the  increased 
thresholds  of  curves  2 and  12  as  compared  to  curve  8 (Figure  49). 

In  17-4  PH  H 1050  a large  difference  (30  ksi  /Tn\  33  MPa  /”m) 
threshold  of  curves  2 and  6 is  seen  if  wave  form  is  considered  con- 
stant and  R-ratio  (0.05  versus  0.67)  is  compared  (Figure  49).  The 
counterpart  in  15-5  PH  but  under  a H 1100  heat  treat  shows  a smaller 
threshold  separation  (6  ksi  ./Tn,  6.6  MPa  ^rii)  for  curves  N and  K. 

The  N and  K curves  do  tend  to  diverge  to  fast  crack  growth  rates 
more  rapid  than  comparable  17-4  PH  curves  2 and  6. 

6.3.3  Wave  Form  Effect 

This  portion  of  the  work  evaluates  the  impact  of  load  wave  form 
on  microstructure  by  considering  alloy  performance  over  an  entire 
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stress  intensity  and  fatigue-crack  growth  range.  Figure  49  can  be 
used  to  compare  10  Hz  sine  wave  to  one  minute  hold-time  trapezoidal 
wave  for  curves  1 and  2,  5 and  6,  8 and  12.  Figure  51  can  be  used 
to  compare  H and  N.  All  comparisons  are  under  the  same  heat  treat, 
R-ratio  and  environment.  Observations  are  that  due  to  substantial 
curve  shifts  in  2 and  12,  17-4  PH  under  this  wave  form  is  responding 
to  what  is  hypothesized  to  be  strain-induced  stress  relaxation 
assisted  by  the  microstructural  feature,  delta  ferrite.  The  premise 
is  that  under  a one  minute  hold-time  at  load,  an  increasing  number 
of  slip  systems  are  activated,  elastically  spreading  applied  strain. 
In  conjunction  with  this  a possibility  exists  that  delta  ferrite 
is  acting  like  a strain  absorber  by  producing  a ductile  rupture  in 
lieu  of  cleaving  In  the  microstructure.  Both  alloys  17-4  PH  and 
15-5  PH  are  basically  similar,  except  for  the  delta  ferrite  present 
in  the  17-4  microstructure.  It  is  suspected  that  delta  ferrite  is 
contributing  to  the  variance  in  fatigue-crack  growth  rate  response 
over  the  entire  AK  range  (Figure  49). 

6.3.4  Environmental  Effect 

This  study  culminates  its  purpose  by  exploring  the  wave  form 
hold-time  effect  in  a 3.5%  ASTM  salt  water  environment.  The  objec- 
tive was  to  relate  laboratory  environment  and  testing  of  a specimen 
with  a comparable  field  environment  and  use  a discrete  wave  form 
element,  relative  to  a structure's  actual  load  spectrum.  Figure  52 
illustrates  a composite  plot  for  the  conditions  of  H 1100,  R = 0.05 
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Figure  52  Composite  Log-Linear  Plot  of  Fatigue-Crack  Growth  Rates  of 
17-4  PH  and  15-5  PH  Used  In  the  Evaluation  of  Environment 
Under  Similar  Heat  Treat  Condition,  R-ratlo  and  Wave  Form 
(12  vs.  10  vs.  N,  N vs.  B,  10  vs.  B,  K vs.  L)  Between 
Materials. 
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and  trapezoidal  one  minute  hold-time  for  17-4  PH  in  ambient  air  and 
salt  water  (curves  12  versus  10)  and  15-5  PH  (curves  N versus  B). 
Here  for  rational  interpretation,  considerations  must  be  three  fold: 
where  are  the  curves  entered  due  to  manufacturing  defects;  what  are 
the  magnitudes  of  deviation  from  a 10  Hz  design  and  decision  base- 
line; and  what  are  stress  ratio  effects  (i.e.  R-ratio,  residual 
stress).  A 15-5  PH  curve  for  H 1100  and  R = 0.67  was  generated. 

The  deviation  from  curve  B (H  1100  and  R s 0.05)  was  substantial 
(Figure  52). 

With  the  above  factors  in  mind,  it  is  believed  from  the  shift- 
ing of  these  curves  for  the  different  test  conditions  that  there  is 
no  doubt  microstructure  dictates  the  alloy  performance  under  fa- 
tigue-crack propagation.  Further,  the  slower  growth  rates  of  17-4 
PH  (curve  10)  between  the  stress  intensity  of  24  ksi  /Tn  (26.4  MPa 
✓m)and  90  ksi  v'Tn  (98.9  MPa  ✓ni)  versus  15-5  PH  (curve  B)  is  be- 
lieved to  be  attributed  to  the  microstructure  being  strain-induced 
to  relax  the  stress  into  the  matrix. 


CHAPTER  7 


MICROSTRUCTURAL  METALLOGRAPHY  AND  FRACTOGRAPHY 

By  collectively  considering  alloy  chemistry,  thermo-mechanical 
history,  metallography  and  fractography  an  explanation  may  be  developed 
for  alloy  performance.  The  rudimentary  factors  requisite  for  under- 
standing mi crostructural  control  to  precisely  predict  component  ma- 
terial life  rest  in  comprehending  the  relation  of  energy  balance, 
molecular  arrangement  and  strain  compatibility  to  the  material  fractur- 
ing process.  Proceeding  with  the  supposition  that  these  factors  are 
e4M'..4l,  coupled  with  verification  of  the  hypothesis  of  this  work 
(i.e.  microstructure  does  dictate  alloy  performance,  chapter  6)  attempts 
are  made  to  discriminate  and  relate  "why"  and  "how"  observed  differ- 
ences between  microstructures  and  between  fracture  surface  determine 
alloy  performance.  Here,  a foundation  exists  for  bridging  material 
performance  to  design  requirements  by  way  of  the  presented  fracture 
and  fatigue  mechanics  methodology  (chapter  5).  This  technique  can 
comprehend  and  decipher  the  effect  of  operational  performance  elements. 

7.1  METALLOGRAPHIC  OBSERVATION 

The  17-4  and  15-5  alloys  were  examined  in  both  heat  treats,  H 
1050  and  H 1100  for  all  three  directions  (T,  S,  L).  At  H 1050  the 
finely  dispersed  copper-rich  particles  are  becoming  observable  at  500x. 
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The  effect  of  these  dispersed  particles  Is  In  preventing  dislocation 
motion  which  depends  In  a complex  way  on  the  particle  size  and  shape, 
conformity  of  the  particle's  crystal  structure  to  the  matrix,  and  the 
elastic  properties  of  the  particle. 

A crystallographic  alignment  between  particles  and  matrix  at  the 
atomic  bonding  level  can  lead  to  complete  coherency.  In  a coherent 
structure  the  particles  produce  elastic  strains  In  the  matrix  develop- 
ing strain  fields  that  resist  free  dislocation  movement  and  are  a major 
source  of  hardening  In  the  alloy.  In  fatigue,  the  copper-rich  precipi- 
tates or  particles  In  17-4  and  15-5  alloys  limit  the  free  dislocation 
slip-path,  (l.e.  across  the  grain)  restricting  dislocation  movement 
over  small  distances.  Thus,  the  presence  of  these  particles  tend  to 
develop  a finer  overall  slip  behavior  requiring  more  energy  and  In- 
creased number  of  repetitive  load  cycles  to  fracture  the  material. 

For  the  H 1100  heat  treat  ,,ondition  particle  size  and  shape  change 
minutely  with  respect  to  particles  surface  area.  The  major  difference 
between  heat  treat  conditions  Is  recognized  In  the  temper  state  of  the 
martensite.  The  H 1100  condition  tends  to  precipitate  more  carbon  In 
the  matrix,  relaxing  Internal  crystal  lattice  strains.  An  optical  com- 
parison can  be  made  between  Figures  53,  54  and  Figures  59,  60  to  see 
the  additional  precipitated  carbon  darkening  the  latter  micrographs 
and  Increasing  the  definition  of  grain  boundaries.  All  micrographs  lend 
themselves  to  a comparison  between  the  heat  treat  condition  as  follows; 

H 1050  Figures  53  through  and  including  58  versus  H 1100  Figures  59 
through  and  Including  64.  Any  retained  austenite  appears  to  be  submicro- 


Figure  53  Microstructure  of  17-4  PH,  H 1050  Heat  Treat  Condition. 

View  Is  Into  the  plane  parallel  to  T In  specimen  T-L, 
Figure  10.  Tempered  martensitic  matrix  with  elongated 
pools  of  delta  ferrite,  and  finely  dispersed  copper-rich 
precipitate  particles.  (Magnification:  500x,  Etch:  Mixed 
Acids). 


Figure  54  Microstructure  of  15-5  PH,  H 1050  Heat  Treat  Condition.  View 
Is  Into  the  plane  parallel  to  T In  specimen  T-L,  Figure  10. 
Tempered  martensitic  matrix,  similar  carbon  presence  as  In 
Figure  53  but  virtually  no  delta  ferrite.  (Magnification: 
500x,  Etch:  Mixed  Acids). 
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Figure  55  Microstructure  of  17-4  PH,  H 1050  Heat  Treat  Condition.  View 
Is  Into  the  plane  parallel  to  L In  specimen  T-L,  Figure  10. 

The  delta  ferrite  pools  are  substantially  Increased  over  those 
In  Figure  53  within  a martensitic  matrix.  (Magnification: 
500x,  Etch:  Mixed  Acids). 


Figure  56  Microstructure  of  15-5  PH,  H 1050  Heat  Treat  Condition.  View 
Is  Into  the  plane  parallel  to  L In  specimen  T-L,  Figure  10. 
Delta  ferrite  appears  to  be  small  oblated  spherlods  through- 
out the  martensitic  matrix.  (Magnification:  500x,  Etch: 

Mixed  Acids). 
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Figure  57  Microstructure  of  17-4  PH,  H 1050  Heat  Treat  Condition.  View 
Is  Into  the  plane  parallel  to  S In  specimen  T-L,  Figure  10. 
The  delta  ferrite  forms  in  well  defined  stringers  In  a marten 
sltlc  matrix.  {Magnification:  500x,  Etch:  Mixed  Acids). 


Figure  58  Microstructure  of  15-5  PH,  H 1050  Heat  Treat  Condition.  View 
Is  Into  the  plane  parallel  to  S In  specimen  T-L,  Figure  10. 
The  matrix  reflects  a primarily  martensitic  structure. 
(Magnification:  500x,  Etch:  Mixed  Acids). 
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Figure  59  Microstructure  of  17-4  PH,  H 1100  Heat  Treat  Condition.  View 
Is  Into  the  plane  parallel  to  T In  specimen  T-L,  Figure  10. 
The  delta  ferrite  stringers  are  prevalent  throughout  the 
martensitic  matrix.  In  this  heat  treat  condition  the  copper- 
rich  finely  dispersed  precipitates  tend  to  Increase  in  size. 
(Magnification:  500x,  Etch:  Mixed  Acids). 


Figure  60  Microstructure  of  15-5  PH,  H 1100  Heat  Treat  Condition.  View 
Is  Into  the  plane  parallel  to  T In  specimen  T-L,  Figure  10. 

The  matrix  tends  to  Increase  Its  carbon  presence  In  the 
martensitic  form  as  the  grain  boundary  definition  Is  enhanced. 
Delta  ferrite  Is  virtually  absent.  (Magnification:  500x, 

Etch:  Mixed  Acids). 


Figure  61  Microstructure  of  17-4  PH,  H 1100  Heat  Treat  Condition.  View 
Is  Into  the  plane  parallel  to  L In  specimen  T-L,  Figure  10. 
The  delta  ferrite,  similar  to  Figure  55,  tends  to  agglomerate 
Into  extended  pools  within  the  martensitic  matrix.  (Magnifi- 
cation: 500x,  Etch:  Mixed  Acids). 


Figure  62  Microstructure  of  15-5  PH,  H 1100  Heat  Treat  Condition.  View 
Is  Into  the  plane  parallel  to  L In  specimen  T-L,  Figure  10. 

The  matrix  has  Increased  carbon  within  the  martensitic  struc- 
ture. Grain  boundaries  are  more  defined.  (Magnification:  500x, 

Etch:  Mixed  Acids). 


Figure  63  Microstructure  of  17-4  PH,  H 1100  Heat  Treat  Condition.  View 
is  Into  the  plane  parallel  to  S in  specimen  T-L,  Figure  10. 

The  structure  has  well  defined  grain  boundaries  due  to  in- 
creased carbon  precipitation  In  H 1100  versus  H 1050.  The 
martensitic  matrix  is  a complex  lamellar  structure  with  copper- 
rich  particles  appearing  throughout.  The  delta  ferrite 
stringers  are  In  the  rolling  direction.  (Magnification:  500x, 

Etch:  Mixed  Acids). 


Figure  64  Microstructure  of  15-5  PH,  H 1100  Heat  Treat  Condition.  View 
is  Into  the  plane  parallel  to  S In  specimen  T-L,  Figure  10. 
The  structure  has  well  defined  grain  boundaries  due  to  In- 
creased carbon  precipitation  in  H 1100  versus  H 1050.  The 
matrix  tends  to  elongate  the  copper-rich  particles.  Delta 
ferrite  is  virtually  non-existent.  (Magnification:  500x, 

Etch:  Mixed  Acids). 
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scoplc  and  not  Influencing  the  measurable  hardness  of  the  alloys. 

The  most  observable  variance  between  17-4  PH  and  15-5  PH  micro- 
structures for  both  heat  treat  conditions  in  all  directions  is  the 
presence  of  delta  ferrite  stringers  in  17-4  PH  and  not  In  15-5  PH. 
Depending  on  the  thickness,  the  stringers  appear  as  either  a white 
or  light  gray  extended  pool  In  micrographs  of  Figures  53,  55,  57,  59, 

61  and  63. 

At  first  glance,  most  Inclusions  are  thought  to  be  incoherent  with 
the  matrix,  having  mismatch  in  elastic  modulus  and  becoming  active 
stress  raisers.  In  general  terms,  inclusions  do  not  have  to  be  detri- 
mental. For  example,  inclusions,  such  as  those  In  free-cutting  steels 
containing  sulphur  and  phosphorus  for  Increased  machinabllity,  do  not 
deteriorate  fatigue  properties  (Forsyth,  1969).  However,  silicate, 
nitride,  hydride,  and  lead  Inclusions  can  be  detrimental  to  fatigue- 
crack  growth. 

As  previously  mentioned  in  section  6.3.4,  the  effect  of  the 
delta  ferrite  Inclusions  in  17-4  PH  appears  to  be  beneficial  over  a 
limited  stress  Intensity  range.  Their  coherency  with  the  matrix  may 
be  adequate  enough  to  produce  a retardation  in  fatigue-crack  growth 
rates  as  the  crack  moves  through  the  microstructure.  Additionally, 
the  delta  ferrite  is  coherent  within  itself  and  may  transfer  strains 
to  more  matrix  by  virtue  of  Its  surface  area,  consequently  Inducing  a 
stress  relaxation  by  Increased  strain  spreading. 


In  Figure  52  the  15-5  (curve  B)  shows  an  initial  resistance  over 
17-4, (curve  10)  by  Indicating  a higher  threshold  va’ue.  Compared  to 
curve  10,  curve  B accelerates  more  rapidly  to  90  ksl  /TrT  (98.9  MPa  /m) 
then  slows  to  lower  crack  growth  rate  than  curve  10  until  final  frac- 
ture. f'  possible  explanation  for  17-4  (curve  10)  performance  may  be 
In  the  presence  of  delta  ferrite. 

Mono tonic  hydrogen  and/or  chlorine  may  diffuse  Into  a material 
under  mechanical  cycling  (e.g.  fatigue)  along  .train  or  stringer  bound- 
aries. The  delta  ferrite  stringers  having  surface  exposures  typical 
of  Figures  54  and  57  in  a manufactured  structure  and  Figures  55  and  58 
for  this  through  cracked  specimen  may  perform  In  one  of  several  ways. 

A possibility  Is  the  stringer  boundaries  may  act  to  a limited  extent 
as  extra  paths  for  atomic  hydroqen  entering  the  matrix.  This  may  re- 
sult In  an  earlier  cracking  stage,  fracturing  small  volumes  of  metal 
quickly.  However,  when  the  applied  strains  Increase,  strain  spreading 
through  the  ferrite  stringers  may  become  more  effective,  perhaps  due 
to  the  entire  volume  and  surface  area  of  the  stringer  now  being  af- 
fected. In  return  the  delta  ferrite  stringer  tends  to  effect  more  of 
the  matrix  surpassing  the  hydrogen  effect.  The  Initial  effect  of  dif- 
fusing atomic  hydrogen  may  not  le  as  Intense  along  the  large  surface 
area  of  the  ferrite  stringers.  El  the  ' one  or  a combination  of  both 
mechanisms  would  permit  the  stringer  phase  to  fracture  In  a manner  of 
microvolt!  coalescence  or  ductile  rupturing.  When  the  stress  Intensity 
increases  to  a high  enough  point  (e.g.  100  ksl  *Tn,  110  MPa  .m)  the 
martensitic  matrix  Influences  the  remaining  life  in  fatigue-crack 
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growth  behavior, 

7.2  FRACTOGRAPHIC  EVALUATION 

The  study  of  fracture  surfaces  provides  a diagnostic  tool  for 
failure  Investigators  and  those  working  In  fracture  mechanics.  The 
haphazard  appearance  can  provide  a discernible  history  of  the  frac- 
ture plus  reveal  something  of  the  basic  mechanism  of  fatigue-crack 
growth. 

Fine  precipitates,  and  undlssolved  carbides,  silicates,  oxides, 
nitrides  or  other  particles  with  a propensity  toward  cleavage  are  fre- 
quently present  In  large  numbers  In  all  but  the  purest  of  metals.  When 
the  matrix  around  these  particles  plastically  deforms,  the  particles 
often  cleave  or  pullout  of  the  matrix  generating  new  surfaces.  These 
surfaces  nucleate  microvoids  that  can  grow.  The  spherical  shape  of 
the  voids  prohibits  simple  glide  and  forces  complex  slip  to  ' rur 
through  these  plastically  worked  surfaces  of  the  void.  As  olds 
grow  and  finally  impinge  on  each  other  the  last  material  separates 
leaving  the  typical  cup-cone  equlaxed  dimple.  If  these  dimples 
coalesce  In  a shearing  direction  they  are  elongated  In  that  direction 
and  called  "shear  dimples".  This  briefly  describes  "microvoid 
coalescence"  or  "dimpled  rupture". 

In  Figure  65  and  Table  X a graphic  representation  is  presented 
for  the  measured  positions  from  which  the  fracture  surfaces  were  sampled. 
The  average  crack  growth  rates  and  stress  Intensities  are  shown  with 
respect  to  these  positions.  All  fracture  specimens  compared  were  tested 


Figure  65  Fractographlc  Layout.  The  average  crack  growth  rate  and 
average  stress  Intensity  are  Indicated  for  the  locations. 
Fracture  sections  were  removed  and  photographed. 
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at  R = 0.05  under  a 10  Hz  sine  wave  In  the  H 1100  condition.  Looking 
at  Specimens  8 and  H (Figures  66  and  67),  a fine  slip  Is  prevalent 
throughout  the  first  two  fracture  sections  with  both  exhibiting  river 
patterns  or  ledges  thought  to  be  signs  of  local  ductile  fracture.  The 
ledges  appear  to  be  at  an  Intermediate  stage  of  growth  when  the  fatigue 
crack  Is  progressing  in  the  "slow  growth"  regime.  The  crack  appears  to 
follow  a jogged  path  Influenced  by  the  local  microstructure  while  break- 
ing In  a ductile  manner  (l.e.  slip)  between  these  various  levels.  In 
the  fast  growth  regime  of  15-5  PH  (Figure  67,  H-3)  the  dimples  are 
finer  and  ledges  tend  to  be  subdued,  indicating  that  in  unstable  frac- 
ture the  sensitivity  to  microstructure  Is  reduced.  Even  though  micro- 
void coalescence  or  dimpled  rupture  does  occur  In  both  alloys  the  17-4 
PH  (Figure  66,  8-3)  does  have  unstable  dimples  on  several  plateaux  that 
are  divided  by  steps,  ridges  or  channels  indicating  microstructural 
sensitivity  under  fast  fracture.  An  explanation  of  15-5  PH  having  a 
higher  fatigue-crack  growth  rate  than  17-4  PH  under  these  conditions 
Is  believed  found  in  the  above  microstructural  effects  on  fracture. 

The  fracture  surfaces  of  the  alloys  tested  under  a one  minute 
hold-time,  Specimens  12  and  N (Figure  68  and  69),  shows  distinct  ledges, 
rivers  and  troughs.  Both  alloys  have  extensive  secondary  cracking  in 
the  longitudinal  and  transverse  depth  directions.  In  Figure  68,  sec- 
tion 12-2,  Intricate  slip  design  occurs  throughout  the  area,  as  the 
main  crack  path  follows  the  microstructure  Intergranularly.  Through  the 
center  of  the  picture  a longitudinal  sequence  of  dimpled  rupture  occurs 


Figure  66  Fracture  Surface  of  Specimen  8,  17-4  PH,  H 1100  Heat  Treat 

Condition,  Tested  under  10  Hz  Sine  Wave,  R = 0.05  In 
. Ambient  Air.  Location  of  fractographs  to  original  fracture 

I surface  with  corresponding  DaDN  and  AK  are  in  Table  X.  The 

| arrow  on  micron  marker  Indicates  crack  growth  direction. 

8-1)  This  view  shows  a slip  design  moving  Intergranularly 
with  modest  cracking  In  the  longitudinal  depth  direc- 
tion. (Magnification:  lOOOx). 

8-2)  Slip  formation  Is  recognizable  with  some  mixed  modes 
of  quasi -cleavage  appearing.  The  long  cracks  In  the 
i direction  of  cracking  appear  to  occur  at  the  bottom 

j of  trenches.  (Magnification:  750x). 

' 8-3)  Dimpled  rupture  In  unstable  fracture  regions  with  1n- 

j tense  pullouts  around  what  appear  to  be  secondary 

j particles.  Well  defined  terraces  or  planes  are  seen 

as  the  crack  moves  through  the  microstructure. 
(Magnification:  750x). 


Figure  67  Fracture  Surface  of  Specimen  H,  15-5  PH,  H 1100  Heat  Treat 
Condition,  Tested  under  10  Hz  Sine  Wave,  R * 0.05  In 
Ambient  Air.  Location  of  fractographs  to  original  fracture 
surface  with  corresponding  DaON  and  AK  are  In  Table  X.  The 
arrow  on  micron  maker  Indicates  crack  growth  direction. 

H-l ) The  view  represents  the  typical  Intergranular  slipped 
surfaces  with  extensive  longitudinal  and  transverse 
depth  cracks,  more  than  In  Specimen  8.  (Magnifi- 
cation: 750x). 

H-2)  Looking  Into  the  depth  of  suspect  Intergranular 

cracking,  slip  ridges  can  be  observed  at  many  levels 
with  these  transverse  depth  cracks.  (Magnification: 
750x) . 

H-3)  Final  unstable  fracture  exhibits  more  dimpled  rupture 
with  noticeable  pullouts  In  conjunction  with  more 
secondary  depth  cracking  than  in  Specimen  8.  (Magni- 
fication: 750x). 
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Figure  68  Fracture  Surface  of  Specimen  12,  17-4  PH,  H 1100  Heat  Treat 
Condition,  Tested  under  a One  Minute  Hold-Time  Trapezoidal 
Wave,  R = 0.05  in  Ambient  Air.  Location  of  fractographs 
to  original  fracture  surface  with  corresponding  DaDN  and  AK 
are  In  Table  X.  The  arrow  on  micron  marker  Indicates  crack 
growth  direction. 

12-1)  Cracking  appears  to  be  more  transgranular  slipping 
with  extensive  secondary  cracking  In  longitudinal 
and  transverse  depth  direction.  (Magnification:  750x). 

12-2)  The  retardation  or  suck-back  of  slip  strlatlons  can  be 
seen  as  they  encounter  a cluster  of  voids  or  pulled 
out  stringers,  believed  to  be  delta  ferrite. 
(Magnification:  750x). 

12-3)  Glide  planes  or  wavy  surfaces,  possibly  serpentine 
glide,  are  visible  at  the  center  right.  In  the 
direction  of  crack  growth,  river  patterns  prevail 
with  slip.  (Magnification:  750x). 


Figure  69  Fracture  Surface  of  Specimen  N,  15-5  PH,  H 1100  Heat  Treat 
Condition,  Tested  under  a One  Minute  Hold-Time  Trapezoidal 
Wave,  R * 0.05  In  Ambient  Air.  Location  of  fractographs 
to  original  fracture  surface  with  corresponding  DaDN  and  AK 
are  In  Table  X.  The  arrow  on  micron  marker  Indicates  crack 
growth  direction. 

N-l ) Intergranular  cracking,  predominantly  slip,  with 
secondary  cracking  and  scattered  pullout  exists. 
(Magnification:  750x). 

N-2)  Slipping  through  a region  with  obvious  mixed  mode 
where  cleavage  Is  a vertical  plane,  concentration 
of  voids  or  pullouts  Is  low,  perhaps  due  to  limited 
delta  ferrite.  (Magnification:  750x). 

N-3)  Ductile  rupture  or  microvoid  coalescence  prevail  as 
dimples  Indicative  of  unstable  ductile  fracture  form 
on  this  area  of  the  specimen.  (Magnification:  750x). 
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as  the  strlatlons  are  sloped  back  or  retarded  with  respect  to  the  crack 
direction.  This  Is  believed  to  be  a delta  ferrite  stringer  rupturing. 

A combination  of  strain  spreading  In  conjunction  with  the  observed 
longitudinal  dimpled  rupture  In  17-4  PH  Is  believed  responsible  for  the 
overall  reduction  of  the  fatigue-crack  growth  rate  and  Increase  In  the 
defined  threshold  as  compared  to  15-5  PH.  Further  15-5  PH  (Figure 
69,  N-2)  exhibits  mixed  mode  cleavage  planes  vertical  to  the  slip  de- 
sign and  crack  direction. 

The  last  sequence  of  fracture  surfaces  represent  the  two  alloys: 
17-4  PH  and  15-5  PH,  (Figures  70  and  71,  respectively)  tested  under 
one  minute  hold-time  In  3.5%  salt  water.  The  17-4  PH  showed  an  Intense 
slip  design  throughout  the  sections.  However,  evident  In  10-1  (Figure 
70)  are  vertical  planes  or  ridges  rising  out  of  the  fracture  surface 
which  may  be  evidence  of  a mechanism  called  slip-plane  initiation, 
persisting  as  part  of  the  growth  process.  This  mechanism  labelled 
Stage  I by  Forsyth  (1969)  is  quickly  succeeded  by  normal  or  Stage  ET 
growth  forming  strlatlons.  Briefly  stated  the  origin  of  fatigue 
cracks  develop  from  slip  bands  extending  along  slip  planes.  The  crack 
origin  Is  associated  with  crystal  slip  or  shear  making  planes  oriented 
approximately  45°  to  the  maximum  tensile  stress  susceptible  to  a maxi- 
mum shear  stress.  There  are  numerous  exceptions  to  this  Idea  but  this 
mode  of  failure  exists.  Cracking  of  this  type  may  be  distinguished  by 
virtually  featureless  fracture  surfaces  or  surface  ridges  In  the  direc- 
tion of  crack  growth  (Forsyth,  1969).  The  remaining  fracture  sections 


Fracture  Surface  of  Specimen  10,  17  4 PH,  H 1100  Heat  Treat 
Condition,  Tested  under  a One  Minute  Hold-Time  Trapezoidal 
jJave,  R = 0.05,  In  3.5t  Salt  Water.  Location  of  fractographs 
to  original  fracture  surfaces  with  corresponding  DaDN  and  &K 
are  *n  Table  X.  The  arrow  on  the  micron  marker  Indicates 
crack  growth  direction. 

10-1)  Majority  of  surface  has  slipped  with  occasional  mixed 
mode  quasi-cleavage  area.  Depth  crack  runs  longitudinal 
to  crack  growth  direction.  (Magnification:  750x). 

10-2)  High  intensity  of  slip,  very  little  mixed  mode.  Intense 
transverse  depth  cracks  and  pullouts  are  seen.  Pullouts 
cr  voids  tend  to  agglomerate  In  a cluster  or  stringer 
formation  a^d  suggest  delta  ferrite  tearing  apart  In  a 
ductile  manner.  (Magnification:  750x). 

10-3)  Moving  to  the  edge  of  unstable  fracture  Initiation, 

slip  design  if  observed  In  the  depth  of  a crater,  mo.,1, 
probably  an  Intergranular,  while  accompanied  by 
secondary  cracking  and  some  minor  mixed  mode.  (Magnlf*- 
catlon:  ’50x). 
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Figure  71  Fracture  Surface  of  Specimen  8,  15-5  PH,  H 1100  Heat  Treat 
Condition,  Tested  under  a One  Minute  Hold-Time  Trapezoidal 
Wave,  R • 0.05,  In  3.5%  Salt  Water.  Location  of  fractographs 
to  original  fracture  surfaces  with  corresponding  OaDN  and  AK 
are  In  Table  X.  The  arrow  on  the  micron  marker  Indicates 
crack  growth  direction. 

6-1)  The  cracking  exhibits  transgranular  cleavage  with 
associate  longitudinal  and  transverse  depth  crack 
In  a mixed  field  of  slip.  (Magnification:  750x). 

B-2)  A typical  cleavage  burst  surrounded  by  quasi - 
cleavage  and  slip  design  at  various  depth  levels. 

The  center  and  the  outside  edges  of  the  burst  appear 
to  be  pullouts  or  voided  points  from  which  the  burst 
may  have  occurred.  (Magnification:  750x). 

B-3)  Mixed  mode  fracture  Is  near-vertical  cleaving  of  the 
microstructure,  then  slipping  and  horizontal  cleaving. 
Notice  pullouts  or  voids  on  the  rounded  slope  In  con- 
junction with  well  defined  strlatlon  design.  (Magnifi- 
cation: 750x). 
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of  Figure  70,  10-2  and  10-3,  show  slipping  In  conjunction  with  sequen- 
tial longitudinal  dimples  rupture.  Exploring  an  Intergranular  bound- 
ary In  10-3,  slip  strlatlons  can  be  observed  with  secondary  cracking 
and  noticeable  pullout  or  void  coalescence.  Specimen  B,  the  15-5  PH 
counterpart  to  17-4  PH,  Specimen  10  clearly  Illustrates  why  Its  fati- 
gue-crack growth  mode  may  be  treacherous  to  a design.  The  material 
with  delta  ferrite  absent  Is  affected  by  the  3.5%  salt  water  environ- 
ment. At  an  early  stage  transgranular  cleavage  fracture  Is  evident 
as  the  applied  stress  Intensity  Increases  with  crack  extension,  cleav- 
age bursting  occurs  from  a central  nucleus.  This  nucleus  or  Initia- 
tion site  for  this  type  of  fracture  can  be  due  to  a cleaved  particle 
resulting  In  a void  or  any  one  of  a number  of  void  generators;  even  a 
"fish  eye"  characteristic  of  hydrogen  embrittlement.  In  the  last 
section  B-3,  vertical  cleavage,  a slipped  or  striated  plateau  and 
another  cleaved  wall  Is  Indicative  of  a mechanism  that  deteriorates 
grain  boundary  cohesion.  By  comparing  the  collection  of  fractures  It 
Is  concluded  that  15-5  PH  Is  more  affected  by  salt  water  and  hold-time 
than  17-4  PH.  Also,  when  the  crack  growth  regimes  enter  a fast  growth 
mode  15-5  PH  relies  on  Its  Inherent  fracture  toughness  to  eluslvely 
suggest  a benefit.  The  transition  point  to  final  fracture  being  an 
area  where  a structures  useful  life  is  almost  ended. 
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CHAPTER  8 • 

CONCLUSIONS  AND  SUMMARY 

8.1  CONCLUSIONS 

As  concluded  from  the  results  of  this  work  and  the  data  herein 
the  following  observations  are  presented. 

. From  this  study  in  accordance  with  the  proposed 
hypothesis  it  is  determined  that  mi crostructural 
control  can  be  used  to  influence  an  alloy's  re- 
sponse in  fatigue-crack  growth. 

. The  results  determined  that  slight  mi crostructural 
alterations  developed  through  heat  treatments 
responded  differently  for  fatigue-crack  growth 
testing  conditions  such  as  R-ratio,  wave  form 
and  environment. 

. A combination  of  salt  water,  hold-time  and  residual 
stress  caused  the  microstructures  to  respond 
differently  than  a 10  Hz  sine  wave,  ambient  air 
fatigue-crack  growth  test. 

. This  study  has  accomplished  the  identification  of 
a unique  hold-time  behavior  in  17-4  PH  not  seen 
nor  even  suspect  before  in  fatigue-crack  growth 
response  of  a martensitic  steel  or  stainless 
steel . 

. A contemporary  Wei  bull -type  mathematical  method 
has  been  demonstrated  to  be  applicable  to  the 
quantification  of  mi crostructural  effects  on  fa- 
tigue-crack growth. 
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. The  data  are  representable  through  four  para- 
meters of  a functional  curve,  not  a straight 
line,  that  precisely  applies  mathematical  boundary 
conditions  to  the  physics  of  the  fatigue-crack 
growth  process. 

Detailed  methods  for  application  of  mathematical  physics  to  fatigue 
crack  growth  are: 


. Exponential  crack  propagation  relations  have  been 
introduced  that  meet  three  key  criteria  of 
acceptability. 

1.  As  AK  approaches  an  upper  limit,  1^,  the  crack 
growth  rate  increases  asymptotically. 

2.  As  AK  approaches  a lower  threshold,  K..,  the 
crack  growth  rate  diminishes  to  zero. 

3.  Slow  and  fast  crack  growth  can  be  distinguished. 

. Mathematical  derived  values  of  stress  intensity 
threshold,  Kth,  can  be  used  to  predict  the  size  of 

the  largest  flaw  which  will  not  grow  for  a given 
loading  condition. 

. The  threshold  stress  intensity,  Kth*  determined  by  e 

and  a calculated  characteristic  value,  v,  the  upper 
limit  of  slow  crack  growth,  can  be  used  as  key 
analytical  initial  conditions  in  planning  in-service 
inspection  schedules. 

. The  included  graphical  representations  of  fatigue 
and  fracture  mechanics  diagrams  provide  a means  to 
consider  load  shedding,  maintenance  and  finite 
life  design  for  field  operational  applications. 
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As  a result  of  this  work,  it  is  believed  from  the  range  of  variables 
covered  by  the  ten  thousand  data  points  generated  herein,  that 
traditional  fatigue-crack  growth  testing  when  combined  with  con- 
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temporary  Weibull  type  fitting  relations  can  be  used  to  predict  com- 
ponent life. 


8.2  SUMMARY 


The  background  work  for  the  study  included  preparation  of 
three  interactive  Conversational  Programming  System  routines,  as 
follows: 


Routine 

POLFIT 


Applicati  }■, 


Conventional  compact  specimen  compliance 
calibration  data  reduction 


WOLFIT  and 
WYBULL 


Derivation  of  Type-A,  -B  or  -C  exponen- 
tial curve  fitting  relations  for  either 
, .iance  calibration  or  fatigue-crack 
propagation  data  reduction 


LIFE 


Prediction  of  crack  length  as  a function 
of  number  of  load  cycles  for  prescribed 
load  history. 


The  utility  of  the  routines  has  been  demonstrated  throughout  chapter 
5 with  the  aid  of  examples  based  on  mill  annealed  Ti-6Al-6V-2Sn, 
high  oxygen  alloy  fracture  mechanics  test  results. 

The  study  technique  was  developed  at  the  Lockheed  Fatigue  and 
Fracture  Mechanics  Laboratory,  Rye  Canyon  Advanced  Design  and  Re- 
search Laboratories  during  the  period  1969-1976  and  has  been  ex- 
panded at  the  University  of  Missouri-Columbia  from  1974  to  the  cur- 
rent data.  This  study  is  the  most  extensive  use  of  the  technique 
specifically  with  regard  to  discriminating  mi crostructural  effects 
through  physically  relevant  curves  in  lieu  of  straight  lines.  The 
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combined  support  of  metallography  and  fractography  In  fatigue-crack 
propagation  testing  has  strengthened  the  analytical  methods  rele- 
vancy and  Interpretation.  Steel  producers  and  specification  writers 
have  empirically  manipulated  mlcrostructural  chemistries  of  these  two 
alloys  to  Increase  weldability  of  17-4  PH  with  the  presence  of  delta 
ferrite,  then  to  Increase  short  transverse  toughness  of  15-5  PH  by 
limiting  delta  ferrite.  These  present  resultant  alloys  when  tested 
for  fatigue-crack  growth  rates  under  variables  that  relate  to  a ve- 
hicle’s operation  profile  can  distinctly  alter  fatigue-crack  growth 
behavior.  Increases  in  fatigue-crack  growth  rates  were  experienced 
with  Increasing  R-ratlo  from  R * 0.05  to  R = 0.67.  This  trend  pre- 
vailed throughout  the  wave  forms  and  heat  treatments  considered. 

The  wave  form  Influenced  a crack  growth  response  substantially 
by  retarding  the  crack  growth  through  a proposed  mlcrostructural 
mechanism  of  strain-induced  stress  relaxation.  Adding  the  salt 
water  environment  to  this  wave  form  effect  reduced  the  threshold  on 
17-4  PH  but  not  on  15-5  PH.  However,  the  fatigue-crack  growth  rates 
on  15-5  PH  accelerated  more  rapidly  than  the  17-4  PH  until  a stress 
Intensity  of  approximately  90  ksl  /Tn  (98.9  MPa  /m)  was  developed. 

An  additional  factor  worthy  of  note  surfaced  during  metallo- 
graphlc  and  fractographlc  analysis.  This  Is,  In  alloys  (e.g.  17-4 
PH  and  15-5  PH)  that  derive  strength  from  precipitated  particles 
blocking  dislocation  motion  during  fatigue  any  application  of  heat 
after  fatiguing  will  cause  increased  precipitation  and  diffusion  In 
the  fatigued  material  due  to  the  stored  energy  from  previous  dlslo- 
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cation  motion  (Forsyth,  1969).  Due  to  these  mechanisms,  if  a struc- 
ture is  cracked,  it  is  believed  to  be  difficult,  if  not  impossible, 
to  weld  repair  for  satisfactory  long  term  operation.  The  heat  in  re- 
pair welding  will  increase  the  size  of  the  particles  for  a zone  of 
material  leading  to  an  incoherent  particle  size  permitting  more  voids 
to  crack  within  the  microstructure.  If  the  structural  component  could 
be  restrained  for  complete  reheat  treatment,  satisfactory  repair  may 
be  possible. 

In  summarizing,  compliance  calibration  and  fatigue-crack  propa- 
gation variables  are  expressible  in  terms  of  fitting  elements  of  a 
distinct  class  of  exponential  curve  fitting  relations. 

Calculation  of  threshold  stress  intensity  range  values,  distinc- 
tion between  slow  and  fast  growth  by  means  of  risk  function  analysis, 
and  life  prediction  by  integration  of  crack  propagation  curve  fitting 
relations  have  been  demonstrated  in  chapter  5.  Except  for  the  last 
integration  all  considered  analytic  methods  have  been  applied  to  the 
data  per  the  test  matrix  In  Table  IX.  Four  parameters  k,  e,  v and 
have  been  used  to  describe  the  test  data  and  microstructural  responses. 

The  microstructure  controls  final  operational  material  response. 

If  a material  does  not  experience  a simulated  or  actual  operational 
profile  of  a structure,  the  resulting  test  data  can  not  be  employed 
to  describe  the  operational  life  of  the  structure. 

Microstructural  sensitivity  exists  in  advanced  high  performance 
structural  materials  and  varies  in  degree  depending  on  the  opera- 
tional conditions  of  interest.  This  uncertainty  of  materials'  re- 
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sponse  in  fatigue-crack  growth  behavior  dictates  a clear  need  for  a 
reality  that  reflects  simulated  or  operational  profiles  in  testing 
and  in  data  representation.  Comprehending  the  need  for  actual  or 

testing  under  the  structure's  conditions  of  Interest  will  culminate 

\ 

in  a worthwhile  and  satisfactory  result  of  structural  integrity.  Ig- 
norance of  this  need  will  end  in  devastation.  There  is  a permanence 
to  this  type  of  consequence  that  ignored  the  responsibility  in  science 
and  engineering.  This  can  best  be  expressed  by  reflecting  on 
Voltaire's  words:  "To  the  living,  one  owes  consideration;  to  the 
dead,  only  the  truth". 
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Fatigue-Crack  Growth  Reduced  Data 


17-4  PH 
H 1050 
R = 0.67 
f = 10  Hz. 
Sine  Wave 
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17-4  PH 
H 1100 
R = 0.05 

One  Minute  Hold-Time 
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